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Abstract 

Over the past three decades, black holes have played an important role in quantum gravity, 
mathematical physics, numerical relativity and gravitational wave phenomenology. However, 
conceptual settings and mathematical models used to discuss them have varied considerably 
from one area to another. Over the last five years a new, quasi-local framework was introduced 
to analyze diverse facets of black holes in a unified manner. In this framework, evolving black 
holes are modelled by dynamical horizons and black holes in equilibrium by isolated horizons. 
We review basic properties of these horizons and summarize applications to mathematical 
physics, numerical relativity, and quantum gravity. This paradigm has led to significant gen- 
eralizations of several results in black hole physics. Specifically, it has introduced a more 
physical setting for black hole thermodynamics and for black hole entropy calculations in 
quantum gravity, suggested a phenomenological model for hairy black holes, provided novel 
techniques to extract physics from numerical simulations, and led to new laws governing the 
dynamics of black holes in exact general relativity. 
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1 Introduction 



Research inspired by black holes has dominated several areas of gravitational physics since the 
early seventies. The mathematical theory turned out to be extraordinarily rich and full of sur- 
prises. Laws of black hole mechanics brought out deep, unsuspected connections between classical 
general relativity, quantum physics, and statistical mechanics |35l I44L 1151 146| . In particular, they 
provided a concrete challenge to quantum gravity which became a driving force for progress in 
that area. On the classical front, black hole uniqueness theorems 16911119) took the community by 
surprise. The subsequent analysis of the detailed properties of Kerr-Newman solutions |65) and 
perturbations thereof (66) constituted a large fraction of research in mathematical general relativ- 
ity in the seventies and eighties. Just as the community had come to terms with the uniqueness 
results, it was surprised yet again by the discovery of hairy black holes jSHllini- Research in this 
area continues to be an active branch of mathematical physics 18l]. The situation has been similar 
in numerical relativity. Since its inception, much of the research in this area has been driven by 
problems related to black holes, particularly their formation through gravitational collapse |153| . 
the associated critical phenomenon ,67i. .105, . and the dynamics leading to their coalescence (see, 
e.g., [Tl 11261 imni 1581 I^ I138| ). Finally, black holes now play a major role in relativistic astrophysics, 
providing mechanisms to fuel the most powerful engines in the cosmos. They are also among the 
most promising sources of gravitational waves, leading to new avenues to confront theory with 
experiments [75] . 

Thus there has been truly remarkable progress on many different fronts. Yet, as the subject 
matured, it became apparent that the basic theoretical framework has certain undesirable features 
from both conceptual and practical viewpoints. Nagging questions have persisted, suggesting the 
need of a new paradigm. 

Dynamical situations 

For fully dynamical black holes, apart from the 'topological censorship' results which restrict 
the horizon topology |ll()l I9()| . there has essentially been only one major result in exact 
general relativity. This is the celebrated area theorem proved by Hawking in the early 
seventies |llll 111^^ : If matter satisfies the null energy condition, the area of the black hole 
event horizon can never decrease. This theorem has been extremely influential because of 
its similarity with the second law of thermodynamics. However, it is a qualitative result; it 
does not provide an explicit formula for the amount by which the area increases in physical 
processes. Now, for a black hole of mass M, angular momentum J, area a, surface gravity 
K, and angular velocity the first law of black hole mechanics, 

6M = -^Sa + nSJ, (1) 

does relate the change in the horizon area to that in the energy and angular momentum, 
as the black hole makes a transition from one equilibrium state to a nearby one |35l 1182) . 
This suggests that there may well be a fully dynamical version of Equation |^ which relates 
the change in the black hole area to the energy and angular momentum it absorbs in fully 
dynamical processes in which the black hole makes a transition from a given state to one 
which is far removed. Indeed, without such a formula, one would have no quantitative control 
on how black holes grow in exact general relativity. Note however that the event horizons 
can form and grow even in a flat region of space-time (see Figure 0] and Section 12.2.21 for 
illustrations). During this phase, the area grows in spite of the fact that there is no flux of 
energy or angular momentum across the event horizon. Hence, in the standard framework 
where the surface of the black hole is represented by an event horizon, it is impossible to 
obtain the desired formula. Is there then a more appropriate notion that can replace event 
horizons? 
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Equilibrium situations 

The zeroth and first laws of black hole mechanics refer to equilibrium situations and small 
departures therefrom. Therefore, in this context, it is natural to focus on isolated black 
holes. It was customary to represent them by stationary solutions of field equations, i.e, 
solutions which admit a time-translational Killing vector field everywhere, not just in a small 
neighborhood of the black hole. While this simple idealization was natural as a starting point, 
it is overly restrictive. Physically, it should be sufficient to impose boundary conditions at 
the horizon which ensure only that the black hole itself is isolated. That is, it should suffice 
to demand only that the intrinsic geometry of the horizon be time independent, whereas the 
geometry outside may be dynamical and admit gravitational and other radiation. Indeed, 
we adopt a similar viewpoint in ordinary thermodynamics; while studying systems such as 
a classical gas in a box, one usually assumes that only the system under consideration is in 
equilibrium, not the whole world. In realistic situations, one is typically interested in the final 
stages of collapse where the black hole has formed and 'settled down' or in situations in which 
an already formed black hole is isolated for the duration of the experiment (see Figure^. In 
such cases, there is likely to be gravitational radiation and non-stationary matter far away 
from the black hole. Thus, from a physical perspective, a framework which demands global 
stationarity is too restrictive. 




Figure 1: Left panel: A typical gravitational collapse. The portion A of the event horizon at late 
times is isolated. Physically, one would expect the first law to apply to A even though the entire 
space-time is not stationary because of the presence of gravitational radiation in the exterior region. 
Right panel: Space-time diagram of a black hole which is initially in equilibrium, absorbs a finite 
amount of radiation, and again settles down to equilibrium. Portions Ai and A2 of the horizon 
are isolated. One would expect the first law to hold on both portions although the space-time is not 
stationary. 

Even if one were to ignore these conceptual considerations and focus just on results, the 
framework has certain unsatisfactory features. Consider the central result, the first law of 
Equation Q. Here, the angular momentum J and the mass M are defined at infinity while 
the angular velocity f2 and surface gravity n are defined at the horizon. Because one has 
to go back and forth between the horizon and infinity, the physical meaning of the first law 
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is not transparent . For instance, there may be matter rings around the black hole which 
contribute to the angular momentum and mass at infinity. Why is this contribution relevant 
to the first law of black hole mechanics? Shouldn't only the angular momentum and mass 
of the black hole feature in the first law? Thus, one is led to ask: Is there a more suitable 
paradigm which can replace frameworks based on event horizons in stationary space-times? 

Entropy calculations 

The first and the second laws suggest that one should assign to a black hole an entropy which 
is proportional to its area. This poses a concrete challenge to candidate theories of quantum 
gravity: Account for this entropy from fundamental, statistical mechanical considerations. 
String theory has had a remarkable success in meeting this challenge in detail for a subclass 
of extremal, stationary black holes whose charge equals mass (the so-called BPS states) |12(J) . 
However, for realistic black holes the charge to mass ratio is less than 10"^^. It has not been 
possible to extend the detailed calculation to realistic cases where charge is negligible and 
matter rings may distort the black hole horizon. From a mathematical physics perspective, 
the entropy calculation should also encompass hairy black holes whose equilibrium states 
cannot be characterized just by specifying the mass, angular momentum and charges at 
infinity, as well as non-minimal gravitational couplings, in presence of which the entropy is no 
longer a function just of the horizon area. One may therefore ask if other avenues are available. 
A natural strategy is to consider the sector of general relativity containing an isolated black 
hole and carry out its quantization systematically. A pre-requisitc for such a program is the 
availability of a manageable action principle and/or Hamiltonian framework. Unfortunately, 
however, if one attempts to construct these within the classical frameworks traditionally used 
to describe black holes, one runs into two difficulties. First, because the event horizon is such 
a global notion, no action principle is known for the sector of general relativity containing 
geometries which admit an event horizon as an internal boundary. Second, if one restricts 
oneself to globally stationary solutions, the phase space has only a finite number of true 
degrees of freedom and is thus 'too small' to adequately incorporate all quantum fluctuations. 
Thus, again, we are led to ask: Is there a more satisfactory framework which can serve as 
the point of departure for a non-perturbative quantization to address this problem? 

Global nature of event horizons 

The future event horizon is defined as the future boundary of the causal past of future null 
infinity. While this definition neatly encodes the idea that an outside observer can not 'look 
into' a black hole, it is too global for many applications. First, since it refers to null infinity, 
it can not be used in spatially compact space-times. Surely, one should be able to analyze 
black hole dynamics also in these space-times. More importantly, the notion is teleological; 
it lets us speak of a black hole only after we have constructed the entire space-time. Thus, 
for example, an event horizon may well be developing in the room you are now sitting in 
anticipation of a gravitational collapse that may occur in this region of our galaxy a million 
years from now. When astrophysicists say that they have discovered a black hole in the center 
of our galaxy, they are referring to something much more concrete and quasi-local than an 
event horizon. Is there a satisfactory notion that captures what they are referring to? 

The teleological nature of event horizons is also an obstruction to extending black hole me- 
chanics in certain physical situations. Consider for example. Figure [3 in which a spherical 
star of mass M undergoes a gravitational collapse. The singularity is hidden inside the null 
surface Ai at r = 2M which is foliated by a family of marginally trapped surfaces and 

^The situation is even more puzzling in the Einstein- Yang-Mills theory where the right side of Equation Q 
acquires an additional term, V5Q. In treatments based on stationary space-times, not only the Yang-Mills charge 
Q, but also the potential V (the analog of Q and k, is evaluated at infinity 11751 . 
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Figure 2: A spherical star of mass M undergoes collapse. Much later, a spherical shell of mass 
SM falls into the resulting black hole. While Ai and A2 are both isolated horizons, only A2 is part 
of the event horizon. 

would be a part of the event horizon if nothing further happens. Suppose instead, after a 
mihion years, a thin spherical shell of mass 5AI collapses. Then Ai would not be a part of 
the event horizon which would actually lie slightly outside Ai and coincide with the surface 
r — 2{M + 6M) in the distant future. On physical grounds, it seems unreasonable to exclude 
Ai a priori from thermodynamical considerations. Surely one should be able to establish the 
standard laws of mechanics not only for the equilibrium portion of the event horizon but also 
for Ai. 

Next, let us consider numerical simulations of binary black holes. Here the main task is to 
construct the space-time containing evolving black holes. Thus, one needs to identify initial 
data containing black holes without the knowledge of the entire space-time and evolve them 
step by step. The notion of a event horizon is clearly inadequate for this. One uses instead 
the notion of apparent horizons (see Section |2.2|) . One may then ask: Can we use apparent 
horizons instead of event horizons in other contexts as well? Unfortunately, it has not been 
possible to derive the laws of black hole mechanics using apparent horizons. Furthermore, as 
discussed in section 2, while apparent horizons are 'local in time' they are still global notions, 
tied too rigidly to the choice of a space-like 3-surface to be directly useful in all contexts. Is 
there a truly quasi-local notion which can be useful in all these contexts? 

Disparate paradigms 

In different communities within gravitational physics, the intended meaning of the term 'black 
hole' varies quite considerably. Thus, in a string theory seminar, the term 'fundamental black 
holes' without further qualification generally refers to the BPS states referred to above - a 
sub-class of stationary, extremal black holes. In a mathematical physics talk on black holes, 
the fundamental objects of interest are stationary solutions to, say, the Einstein-Higgs- Yang- 
Mills equations for which the uniqueness theorem fails. The focus is on the ramifications 
of 'hair', which are completely ignored in string theory. In a numerical relativity lecture, 
both these classes of objects are considered to be so exotic that they are excluded from 
discussion without comment. The focus is primarily on the dynamics of apparent horizons 
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in general relativity. In astrophysically interesting situations, the distortion of black holes by 
external matter rings, magnetic fields and other black holes is often non- negligible [8711981 155] , 
While these illustrative notions seem so different, clearly there is a common conceptual core. 
Laws of black hole mechanics and the statistical mechanical derivation of entropy should go 
through for all black holes in equilibrium. Laws dictating the dynamics of apparent horizons 
should predict that the final equilibrium states are those represented by the stable stationary 
solutions of the theory. Is there a paradigm that can serve as an unified framework to 
establish such results in all these disparate situations? 

These considerations led to the development of a new, quasi-local paradigm to describe black 
holes. This framework was inspired by certain seminal ideas introduced by Hayward [1171 11181 
11151 1116) in the mid-nineties and has been systematically developed over the past five years or so. 
Evolving black holes are modelled by dynamical horizons while those in equilibrium are modelled 
by isolated horizons. Both notions are quasi-local. In contrast to event horizons, neither notion 
requires the knowledge of space-time as a whole or refers to asymptotic fiatness. Furthermore, 
they are space-time notions. Therefore, in contrast to apparent horizons, they are not tied to 
the choice of a partial Cauchy slice. This framework provides a new perspective encompassing all 
areas in which black holes feature: quantum gravity, mathematical physics, numerical relativity, 
and gravitational wave phenomenology. Thus, it brings out the underlying unity of the subject. 
More importantly, it has overcome some of the limitations of the older frameworks and also led to 
new results of direct physical interest. 

The purpose of this article is to review these developments. The subject is still evolving. Many 
of the key issues are still open and new results are likely to emerge in the coming years. Nonetheless, 
as the Editors pointed out, there is now a core of results of general interest and, thanks to the 
innovative style of Living Reviews, we will be able to incorporate new results through periodic 
updates. 

Applications of the quasi- local framework can be summarized as follows: 

Black hole mechanics 

Isolated horizons extract from the notion of Killing horizons, just those conditions which 
ensure that the horizon geometry is time independent; there may be matter and radiation 
even nearby [HHI- Yet, it has been possible to extend the zeroth and first laws of black hole 
mechanics to isolated horizons 26. .75. .15, . Furthermore, this derivation brings out a con- 
ceptually important fact about the first law. Recall that, in presence of internal boundaries, 
time evolution need not be Hamiltonian (i.e., need not preserve the symplectic structure). If 
the inner boundary is an isolated horizon, a necessary and sufficient condition for evolution 
to be Hamiltonian turns out to be precisely the first law! Finally, while the first law has the 
same form as before (Equation 1^), all quantities which enter the statement of the law now 
refer to the horizon itself. This is the case even when non-Abelian gauge fields are included. 

Dynamical horizons allow for the horizon geometry to be time dependent. This framework 
has led to a quantitative relation between the growth of the horizon area and the flux of energy 
and angular momentum across it |3UI 13 1| . The processes can be in the non-linear regime of 
exact general relativity, without any approximations. Thus, the second law is generalized 
and the generalization also represents an integral version of the first law , applicable also 
when the black hole makes a transition from one state to another, which may be far removed. 

Quantum gravity 

The entropy problem refers to equilibrium situations. The isolated horizon framework pro- 
vides an action principle and a Hamiltonian theory which serves as a stepping stone to 
non-perturbative quantization. Using the quantum geometry framework, a detailed theory 
of the quantum horizon geometry has been developed. The horizon states are then counted 
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to show that the statistical mechanical black hole entropy is indeed proportional to the 
area fTUlimiMI 11491 05] , This derivation is applicable to ordinary, astrophysical black holes 
which may be distorted and far from extremality. It also encompasses cosmological horizons 
to which thermodynamical considerations are known to apply |99| . Finally, the arena for this 
derivation is the curved black hole geometry, rather than a system in flat space-time which 
has the same number of states as the black hole |174l ll4iT] . Therefore, this approach has a 
greater potential for analyzing physical processes associated with the black hole. 

The dynamical horizon framework has raised some intriguing questions about the relation 
between black hole mechanics and thermodynamics in fully dynamical situations In 
particular, they provide seeds for further investigations of the notion of entropy in non- 
equilibrium situations. 

Mathematical physics 

The isolated horizon framework has led to a phenomenological model to understand properties 
of hairy black holes |5T1I2D|- In this model, the hairy black hole can be regarded as a bound 
state of an ordinary black hole and a soliton. A large number of facts about hairy black holes 
had accumulated through semi-analytical and numerical studies. Their qualitative features 
are explained by the model. 

The dynamical horizon framework also provides the groundwork for a new approach to Pen- 
rose inequalities which relate the area of cross-sections of the event horizon Ac on a Cauchy 
surface with the ADM mass Madm at infinity |157| : Ac/16tt < Madm- Relatively recently, 
the conjecture has been proved in time symmetric situations. The basic monotonicity formula 
of the dynamical horizon framework could provide a new avenue to extend the current proofs 
to non-time-symmetric situations. It may also lead to a stronger version of the conjecture 
where the ADM mass is replaced by the Bondi mass 

Numerical relativity 

The framework has provided a number of tools to extract physics from numerical simulations 
in the near-horizon, strong field regime. First, there exist expressions for mass and angu- 
lar momentum of dynamical and isolated horizons which enable one to monitor dynamical 
processes occurring in the simulations and extract properties of the final equilibrium 
state ^IISS]. These quantities can be calculated knowing only the horizon geometry and do 
not pre-suppose that the equilibrium state is a Kerr horizon. The computational resources 
required in these calculations are comparable to those employed by simulations using cruder 
techniques, but the results are now invariant and interpretation is free from ambiguities. Re- 
cent work |34| has shown that these methods are also numerically more accurate and robust 
than older ones. 

Surprisingly, there are simple local criteria to decide whether the geometry of an isolated 
horizon is that of the Kerr horizon I 142j . These criteria have already been implemented 
in numerical simulations. The isolated horizon framework also provides invariant, practical 
criteria to compare near-horizon geometries of different simulations (12) and leads to a new 
approach to the problem of extracting wave-forms in a gauge invariant fashion. Finally, the 
framework provides natural boundary conditions for the initial value problem for black holes 
in quasi-equilibrium |72l 11241 181| . and to interpret certain initial data sets (136) . Many of 
these ideas have already been implemented in some binary black hole codes [SSI |^ 02] and 
the process is continuing. 

Gravitational wave phenomenology 

The isolated horizon framework has led to a notion of horizon multipole moments j24j . They 
provide a diffeomorphism invariant characterization of the isolated horizon geometry. They 
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are distinct from the Hansen multipoles in stationary space-times |1()7) normally used in the 
analysis of equations of motion because they depend only on the isolated horizon geometry 
and do not require global stationarity. They represent source multipoles rather than Hansen's 
field multipoles. In Kerr space-time, while the mass and angular momenta agree in the two 
regimes, quadrupole moments do not; the difference becomes significant when a ~ M, i.e., in 
the fully relativistic regime. In much of the literature on equations of motion of black holes, 
the distinction is glossed over largely because only field multipoles have been available in the 
literature. However, in applications to equations of motion, it is the source multipoles that 
are more relevant, whence the isolated horizon multipoles are likely to play a significant role. 

The dynamical horizon framework enables one to calculate mass and angular momentum of 
the black hole as it evolves. In particular, one can now ask if the black hole can be first 
formed violating the Kerr bound a < M but then eventually settle down in the Kerr regime. 
Preliminary considerations fail to rule out this possibility, although the issue is still open '31J. 
The issue can be explored both numerically and analytically. The possibility that the bound 
can indeed be violated initially has interesting astrophysical implications 

In this review, we will outline the basic ideas underlying dynamical and isolated horizon frame- 
works and summarize their applications listed above. The material is organized as follows. In 
Section 13 we recall the basic definitions, motivate the assumptions and summarize their implica- 
tions. In Section |31 we discuss the area increase theorem for dynamical horizons and show how it 
naturally leads to an expression for the flux of gravitational energy crossing dynamical horizons. 
Section 2] is devoted to the laws of black hole mechanics. We outline the main ideas using both 
isolated and dynamical horizons. In the next three sections we review applications. Section El 
summarizes applications to numerical relativity. Sectional to black holes with hair, and Section d 
to the quantum entropy calculation. Section|Sldiscusses open issues and directions for future work. 
Having read Section\^ Sections\^^\^\^ and^ are fairly self contained and the three applications 
can be read independently of each other. 

All manifolds will be assumed to be C'^"'"^ (with fc > 3) and orientable, the space-time metric 
will be C*^, and matter fields C^^^ . For simplicity we will restrict ourselves to 4-dimensional space- 
time manifolds M (although most of the classical results on isolated horizons have been extended 
to 3-dimensions space-times |^, as well as higher dimensional ones |141p . The space-time metric 
Qab has signature (— , +, -|-, and its derivative operator will be denoted by V. The Riemann 
tensor is defined by Rabc^Wd ■— 2V[aVb]Wc, the Ricci tensor by Rat ■— Racb'^, and the scalar 
curvature by i? := g°'^Rab- We will assume the field equations 

Rab - ^Rgab + ^gab = SnGTab- (2) 

(With these conventions, de Sitter space-time has positive cosmological constant A.) We assume 
that Tab satisfies the dominant energy condition (although, as the reader can easily tell, several 
of the results will hold under weaker restrictions.) Cauchy (and partial Cauchy) surfaces will be 
denoted by M, isolated horizons by A, and dynamical horizons by H . 
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2 Basic Notions 



This section is divided into two parts. The first introduces isolated horizons, and the second 
dynamical horizons. 

2.1 Isolated horizons 

In this part, we provide the basic definitions and discuss geometrical properties of non-expanding, 
weakly isolated, and isolated horizons which describe black holes which are in equilibrium in an 
increasingly stronger sense. 

These horizons model black holes which are themselves in equilibrium, but in possibly dynamical 
space-times j^lEUHHlE]- For early references with similar ideas, see |156I I106'. A useful example 
is provided by the late stage of a gravitational collapse shown in Figure ^ In such physical 
situations, one expects the back-scattered radiation falling into the black hole to become negligible 
at late times so that the 'end portion' of the event horizon (labelled by A in the figure) can be 
regarded as isolated to an excellent approximation. This expectation is borne out in numerical 
simulations where the backscattering effects typically become smaller than the numerical errors 
rather quickly after the formation of the black hole (see, e.g., |34l 148) 1. 

2.1.1 Definitions 

The key idea is to extract from the notion of a Killing horizon the minimal conditions which 
are necessary to define physical quantities such as the mass and angular momentum of the black 
hole and to establish the zeroth and the first laws of black hole mechanics. Like Killing horizons, 
isolated horizons are null, 3-dimensional sub-manifolds of space-time. Let us therefore begin by 
recalling some essential features of such sub-manifolds, which we will denote by A. The intrinsic 
metric qat on A has signature (0,-t-,-|-), and is simply the pull-back of the space-time metric to A, 
Qab — ^ab, whcrc an underarrow indicates the puUback to A. Since qab is degenerate, it does not 
have an inverse in the standard sense. However, it does admit an inverse in a weaker sense: g"'' 
will be said to be an inverse of Qab if it satisfies qamlbnl"^"' = Qab- As one would expect, the inverse 
is not unique: We can always add to q""^ a term of the type £^°-V''\ where is a null normal to A 
and any vector field tangential to A. All our constructions will be insensitive to this ambiguity. 
Given a null normal to A, the expansion is defined as 

e(,) :- q'^'Vah- (3) 

(Throughout this review, we will assume that is future directed.) We can now state the first 
definition: 

Definition 1: A sub-manifold A of a space-time {M,gab) is said to he a non-expanding horizon 
(NEH) if 

1. A is topologically S*^ x M and null; 

2. any null normal of A has vanishing expansion, = 0; and 

3. all equations of motion hold at A and the stress energy tensor Tab is such that —T^^£^ is 
future-causal for any future directed null normal £"• . 

The motivation behind this definition can be summarized as follows. Condition^is imposed for 
definiteness; while most geometric results are insensitive to topology, the 5^ x M case is physically 
the most relevant one. Condition 13 is satisfied by all classical matter fields of direct physical 
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interest. The key condition in the above definition is Condition |21 which is equivalent to requiring 
that every cross-section of A be marginally trapped. (Note incidentally that if 0(f) vanishes for one 
null normal £° to A, it vanishes for all.) Condition|21is equivalent to requiring that the infinitesimal 
area element is Lie dragged by the null normal £"'. In particular, then, Condition |2] implies that 
the horizon area is 'constant in time'. We will denote the area of any cross section of A by oa and 
define the horizon radius as Ra '■= \J aA/47r. 

Because of the Raychaudhuri equation. Condition |21 also implies 

-Rafcr^'' + aabfj"'' = 0, (4) 

where a°-^ is the shear of £a, defined by <Jah '■= ^(a^b) ~ ^^{e)1ab, where the underarrow denotes 
'pull-back to A'. Now the energy condition O implies that Rabi"'^'' is non-negative, whence we 
conclude that each of the two terms in the last equation vanishes. This in turn implies that 
Tabi^ = and W(a£b) = on A. The first of these equations constrains the matter fields on A in 
an interesting way, while the second is equivalent to CiQab = on A. Thus, the intrinsic metric on 
an NEH is 'time-independent'; this is the sense in which an NEH is in equilibrium. 

The zeroth and first laws of black hole mechanics require an additional structure, which is 
provided by the concept of a weakly isolated horizon. To arrive at this concept, let us first 
introduce a derivative operator I? on A. Because qab is degenerate, there is an infinite number 
of (torsion-free) derivative operators which are compatible with it. However, on an NEH, the 
property 'V(aib) = implies that the space-time connection V induces a unique (torsion-free) 
derivative operator 2? on A which is compatible with qab Weakly isolated horizons are 

characterized by the property that, in addition to the metric Qab, the connection component "Da^^ 
is also 'time independent'. 

Two null normals and to an NEH A are said to belong to the same equivalence class [£] 
if = c£" for some positive constant c. Then, weakly isolated horizons are defined as follows: 

Definition 2: The pair (A, [£]) is said to constitute a weakly isolated horizon (WIH) provided 
A is an NEH and each null normal £"" in [£] satisfies 

{CtVa - VaCi)£'' = 0. (5) 



It is easy to verify that every NEH admits null normals satisfying Equation |(SJ, i.e., can be 
made a WIH with a suitable choice of [£]. However the required equivalence class is not unique, 
whence an NEH admits distinct WIH structures 1161. 

Compared to conditions required of a Killing horizon, conditions in this definition are very 
weak. Nonetheless, it turns out that they are strong enough to capture the notion of a black hole 
in equilibrium in applications ranging from black hole mechanics to numerical relativity. (In fact, 
many of the basic notions such as the mass and angular momentum are well-defined already on 
NEHs although intermediate steps in derivations use a WIH structure.) This is quite surprising 
at first because the laws of black hole mechanics were traditionally proved for globally stationary 
black holes |182) . and the definitions of mass and angular momentum of a black hole first used in 
numerical relativity implicitly assumed that the near horizon geometry is isometric to Kerr .5^. 

Although the notion of a WIH is sufficient for most applications, from a geometric viewpoint, 
a stronger notion of isolation is more natural: The full connection V should be time-independent. 
This leads to the notion of an isolated horizon. 

Definition 3: A WIH (A, [£]) is said to constitute an isolated horizon (IH) if 

(CeVa - VaCi)V'' = (6) 
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for arbitrary vector fields tangential to A. 

While an NEH can always be given a WIH structure simply by making a suitable choice of 
the null normal, not every WIH admits an IH structure. Thus, the passage from a WIH to an 
IH is a genuine restriction ^H]- However, even for this stronger notion of isolation, conditions in 
the definition are local to A. Furthermore, the definition only uses quantities intrinsic to A; there 
are no restrictions on components of any fields transverse to A. (Even the full 4- metric gab need 
not be time independent on the horizon.) Robinson-Trautman solutions provide explicit examples 
of isolated horizons which do not admit a stationary Killing field even in an arbitrarily small 
neighborhood of the horizon "HS]. In this sense, the conditions in this definition are also rather 
weak. One expects them to be met to an excellent degree of approximation in a wide variety of 
situations representing late stages of gravitational collapse and black hole mergers^. 

2.1.2 Examples 

The class of space-times admitting NEHs, WIHs, and IHs is quite rich. First, it is trivial to verify 
that any Killing horizon which is topologically S'^ x K is also an isolated horizon. This in particular 
implies that the event horizons of all globally stationary black holes, such as the Kerr-Newman 
solutions (including a possible cosmological constant), are isolated horizons. (For more exotic 
examples, see 11551 .1 



s 




Figure 3: Set-up of the general characteristic initial value formulation. The Weyl tensor component 
^/q on the null surface A is part of the free data which vanishes if A is an IH. 

But there exist other non-trivial examples as well. These arise because the notion is quasi-local, 
referring only to fields defined intrinsically on the horizon. First, let us consider the sub-family of 
Kastor-Traschen solutions |125II152] which are asymptotically de Sitter and admit event horizons. 
They are interpreted as containing multiple charged, dynamical black holes in presence of a positive 
cosmological constant. Since these solutions do not appear to admit any stationary Killing fields, 
no Killing horizons are known to exist. Nonetheless, the event horizons of individual black holes are 
WIHs. However, to our knowledge, no one has checked if they are IHs. A more striking example 
is provided by a sub-family of Robinson-Trautman solutions, analyzed by Chrusciel |68| . These 
space-times admit IHs whose intrinsic geometry is isomorphic to that of the Schwarzschild isolated 
horizons but in which there is radiation arbitrarily close to A. 

■^However, Condition may be too strong in some problems, e.g., in the construction of quasi-equilibrium initial 
data sets, where the notion of WIH is more useful fTM (see Section Ot. 
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More generally, using the characteristic initial value formulation |92l 1161) . Lewandowski jl4()| 
has constructed an infinite dimensional set of local examples. Here, one considers two null surfaces 
A and Af intersecting in a 2-sphere S (see Figure . One can freely specify certain data on these 
two surfaces which then determines a solution to the vacuum Einstein equations in a neighborhood 
of S bounded by A and J\f, in which A is an isolated horizon. 

2.1.3 Geometrical properties 

Rescaling freedom in l"" 

As we remarked in Section f2.1.1l there is a functional rescaling freedom in the choice of a null 
normal on an NEH and, while the choice of null normals is restricted by the weakly isolated 
horizon condition ©, considerable freedom still remains. That is, a given NEH A admits an 
infinite number of WIH structures (A, [£]) 16 . 

On IHs, by contrast, the situation is dramatically different. Given an IH (A, [£]), generically 
the Condition JHl in Definition |31 can not be satisfied by a distinct equivalence class of null 
normals [£']. Thus on a generic IH, the only freedom in the choice of the null normal is that 
of a rescaling by a positive constant 16 . This freedom mimics the properties of a Killing 
horizon since one can also rescale the Killing vector by an arbitrary constant. The triplet 
{qabTT^a, [^°]) is said to constitute the geometry of the isolated horizon. 

Surface gravity 

Let us begin by defining a 1-form uja which will be used repeatedly. First note that, by 
Definition ^ is expansion free and shear free. It is automatically twist free since it is a 
normal to a smooth hypersurface. This means that the contraction of Va^6 with any two 
vectors tangent to A is identically zero, whence there must exist a 1-form uja on A such that 
for any V"" tangent to A, 

VVai'' = V'LUae''. (7) 

Note that the WIH condition lO requires simply that uja be time independent, CiOJa — 0. 
Given Wq, the surface gravity K(i) associated with a null normal is defined as 

Twa. (8) 

Thus, is simply the acceleration of £"'. Note that the surface gravity is not an intrinsic 
property of a WIH (A, [£]). Rather, it is a property of a null normal to A: K(c£) = CK(^(y An 
isolated horizon with k^^) = is said to be an extremal isolated horizon. Note that while the 
value of surface gravity refers to a specific null normal, whether a given WIH is extremal or 
not is insensitive to the permissible rescaling of the normal. 

Curvature tensors on A 

Consider any (space-time) null tetrad (£", n", to", to") on A such that i"" is a null normal to 
A. Then, it follows from Definition ^ that two of the Newman-Penrose Weyl components 
vanish on A: ^-q := Cabcd^m^e'^m'^ = and ^-i := Cabcd'"'m!'(-''n'^ = 0. This in turn implies 
that 'I'2 := Cabcd°''m^'frfn'^ is gauge invariant (i.e., does not depend on the specific choice of 
the null tetrad satisfying the condition stated above.) The imaginary part of is related 
to the curl of uJa, 

duj^2 (Im^'a) e, (9) 

where tab is the natural area 2-form on A. Horizons on which Im5'2 vanishes are said to be 
non-rotating: On these horizons all angular momentum multipoles vanish |24| . Therefore, 
Im4'2 is sometimes referred to as the rotational scalar and as the rotation 1-form of the 
horizon. 
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Next, let us consider the Ricci-tensor components. On any NEH A we have: $00 := 
\Rab^°'(-^ := 0, $01 := ^Rab(^°'m'' = 0. In the Einstein-Maxwell theory, one further has: 
On A, $02 := ^Rabrn'^rJ = and $20 := \Rabm''m!' = 0. 

Free data on an isolated horizon 

Given the geometry (qabj'C), [£]) of an IH, it is natural to ask for the minimum amount of 
information, i.e., the free data, required to construct it. This question has been answered in 
detail (also for WIHs) ^B]. For simplicity, here we will summarize the results only for the 
non-extremal case. (For the extremal case, see 1143) .') 

Let S* be a spherical cross section of A. The degenerate metric qab naturally projects to 
a Riemannian metric qab on S, and similarly the 1-form uja of Equation (|7|) projects to a 
1-form Wo on S. If the vacuum equations hold on A, then given {(jab^'^a) on S, there is, 
up to diffeomorphisms, a unique non-extremal isolated horizon geometry {qab,T^,[i]) such 
that qab is the projection of qab, ^a is the projection of the oja constructed from 2?, and 
= 7^ 0. (If the vacuum equations do not hold, the additional data required is the 
projection on S of the space-time Ricci tensor.) 

The underlying reason behind this result can be sketched as follows. First, since qab is 
degenerate along its non-trivial part is just its projection qab- Second, qab fixes the 
connection on S\ it is only the quantity Sab '■= T^a^b that is not constrained by qab, where 
Ua is a 1-form on A orthogonal to S, normalized so that t°'na = — 1. It is easy to show 
that Sab is symmetric and the contraction of one of its indices with I"' gives uja'- (-""Sab = ^a- 
Furthermore, it turns out that if w^^" ^ 0, the field equations completely determine the 
angular part of Sab in terms of uja and qab- Finally, recall that the surface gravity is not 
fixed on A because of the rescaling freedom in thus the ^-component of LOa is not part of 
the free data. Putting all these facts together, we see that the pair {qab,^a) enables us to 
reconstruct the isolated horizon geometry uniquely up to diffeomorphisms. 

Rest frame of a non-expanding horizon 

As at null infinity, a preferred foliation of A can be thought of as providing a 'rest frame' 
for an isolated horizon. On the Schwarzschild horizon, for example, the 2-spheres on which 
the Eddington-Finkelstein advanced time coordinate is constant - which are also integral 
manifolds of the rotational Killing fields - provide such a rest frame. For the Kerr metric, 
this foliation generalizes naturally. The question is whether a general prescription exists to 
select such a preferred foliation. 

On any non-extremal NEH, the 1-form LOa can be used to construct preferred foliations of A. 
Let us first examine the simpler, non-rotating case in which Im ^2 = 0. Then Equation © 
implies that tOa is curl-free and therefore hypersurface orthogonal. The 2-surfaces orthogonal 
to uJa must be topologically S^ because, on any non-extremal horizon, €°Wa 7^ 0. Thus, 
in the non-rotating case, every isolated horizon comes equipped with a preferred family of 
cross-sections which defines the rest frame |26) . Note that the projection Qa of tOa on any 
leaf of this foliation vanishes identically. 

The rotating case is a little more complicated since cua is then no longer curl-free. Now the 
idea is to exploit the fact that the divergence of the projection uja of uJa on a cross-section is 
sensitive to the choice of the cross-section, and to select a preferred family of cross-sections 
by imposing a suitable condition on this divergence jl6j . A mathematically natural choice 
is to ask that this divergence vanish. However, (in the case when the angular momentum is 
non-zero) this condition does not pick out the v — const, cuts of the Kerr horizon where v is 
the (Carter generalization of the) Eddington-Finkelstein coordinate. Pawlowski has provided 
another condition that also selects a preferred foliation and reduces to the v — const, cuts of 
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the Kerr horizon: 



divw = -Aln |\I'2|^ 



(10) 



where A is the Laplacian of qab- On isolated horizons on which 1^*2 1 is nowhere zero - a 
condition satisfied if the horizon geometry is 'near' that of the Kerr isolated horizon - this 
selects a preferred foliation and hence a rest frame. This construction is potentially useful to 
numerical relativity. 

Symmetries of an isolated horizon 

By definition, a symmetry of an IH (A, [£]) is a diffeomorphism of A which preserves the 
geometry (got, 2^, [^])- (On a WIH, the symmetry has to preserve (gafc, Wa, [^]). There are 
again three universality classes of symmetry groups as on an IH.) Let us denote the symmetry 
group by Ga- First note that diffeomorphisms generated by the null normals in are 
symmetries; this is already built into the very definition of an isolated horizon. The other 
possible symmetries are related to the cross-sections of A. Since we have assumed the cross- 
sections to be topologically spherical and since a metric on a sphere can have either exactly 
three, one or zero Killing vectors, it follows that Ga can be of only three types |15j: 

• Type I: The pair {qab,T^a) is spherically symmetric; Ga is four dimensional. 

• Type II: The pair {qabiT^a) is axisymmetric; Ga is two dimensional. 

• Type III: Diffeomorphisms generated by are the only symmetries of the pair [qab, T^a)] 
Ga is one dimensional. 

In the asymptotically flat context, boundary conditions select a universal symmetry group 
at spatial infinity, e.g., the Poincare group, because the space-time metric approaches a fixed 
Minkowskian one. The situation is completely different in the strong field region near a black 
hole. Because the geometry at the horizon can vary from one space-time to another, the 
symmetry group is not universal. However, the above result shows that the symmetry group 
can be one of only three universality classes. 

2.2 Dynamical horizons 

This section is divided into three parts. In the first, we discuss basic definitions, in the second we 
introduce an explicit example, and in the third we analyze the issue of uniqueness of dynamical 
horizons and their role in numerical relativity. 

2.2.1 Definitions 

To explain the evolution of ideas and provide points of comparison, we will introduce the notion 
of dynamical horizons following a chronological order. Readers who are not familiar with causal 
structures can go directly to Definition[Slof dynamical horizons (for which a more direct motivation 
can be found in [ST]). 

As discussed in Section^ while the notion of an event horizon has proved to be very convenient 
in mathematical relativity, it is too global and teleological to be directly useful in a number 
of physical contexts ranging from quantum gravity to numerical relativity to astrophysics. This 
limitation was recognized early on (see, e.g., page 319) and alternate notions were introduced 

to capture the intuitive idea of a black hole in a quasi-local manner. In particular, to make the 
concept 'local in time'. Hawking |llll 1113] introduced the notions of a trapped region and an 
apparent horizon, both of which are associated to a space-like 3-surface M representing 'an instant 
of time'. Let us begin by recalling these ideas. 
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Hawking's outer trapped surface 5 is a compact, space-like 2-diniensional sub-manifold in 
{A4,gab) such that the expansion 0(£-) of the outgoing null normal to S is non-positive. Hawking 
then defined the trapped region T{M) in a surface M as the set of all points in M through which 
there passes an outer-trapped surface, lying entirely in M. Finally, Hawking's apparent horizon 
dT{M) is the boundary of a connected component of T{M). The idea then was to regard each 
apparent horizon as the instantaneous surface of a black hole. One can calculate the expansion 
G(£) of S knowing only the intrinsic 3-metric qab and the extrinsic curvature Kab of M . Hence, to 
find outer trapped surfaces and apparent horizons on M, one does not need to evolve {qab, Kab) 
away from M even locally. In this sense the notion is local to M . However, this locality is achieved 
at the price of restricting S to lie in M . If we wiggle M even slightly, new outer trapped surfaces 
can appear and older ones may disappear. In this sense, the notion is still very global. Initially, 
it was hoped that the laws of black hole mechanics can be extended to these apparent horizons. 
However, this has not been possible because the notion is so sensitive to the choice of M . 

To improve on this situation, in the early nineties Hayward proposed a novel modification of this 
framework I117f . The main idea is to free these notions from the complicated dependence on M. 
He began with Penrose's notion of a trapped surface. A trapped surface 5* a la Penrose is a compact, 
space-like 2-dimensional sub-manifold of space-time on which 8(£)8(„) > 0, where and are 
the two null normals to S. We will focus on future trapped surfaces on which both expansions are 
negative. Hayward then defined a space-time trapped region. A trapped region T a la Hayward 
is a subset of space-time through each point of which there passes a trapped surface. Finally, 
Hayward's trapping boundary dT is a connected component of the boundary of an inextendible 
trapped region. Under certain assumptions (which appear to be natural intuitively but technically 
are quite strong), he was able to show that the trapping boundary is foliated by marginally trapped 
surfaces (MTSs), i.e., compact, space-like 2-dimensional sub-manifolds on which the expansion of 
one of the null normals, say vanishes and that of the other, say n°, is everywhere non-positive. 
Furthermore, £„8(£) is also everywhere of one sign. These general considerations led him to define 
a quasi-local analog of future event horizons as follows: 

Definition 4: A future, outer, trapping horizon (FOTH) is a smooth 3-dimensional sub-manifold 
H of space-time, foliated by closed 2-manifolds S, such that 

1. the expansion of one future directed null normal to the foliation, say d"" , vanishes, = 0; 

2. the expansion of the other future directed null normal is negative, Q(n) < 0; and 

3. the directional derivative ofQ(^g) along n° is negative, £„ 0(£) < 0. 

In this definition. Condition |2 captures the idea that His a. future horizon (i.e., of black hole 
rather than white hole type), and Condition El encodes the idea that it is 'outer' since infinitesimal 
motions along the 'inward' normal n° makes the 2-surface trapped. (Condition also serves to 
distinguish black hole type horizons from certain cosmological ones |117) which are not ruled out 
by Condition 12) . Using the Raychaudhuri equation, it is easy to show that _ff is either space-like 
or null, being null if and only if the shear aab of as well as the matter flux Tab(°'£'' across H 
vanishes. Thus, when His null, it is a non-expanding horizon introduced in Section f2.1l Intuitively, 
H is space-like in the dynamical region where gravitational radiation and matter fields are pouring 
into it and is null when it has reached equilibrium. 

In truly dynamical situations, then, H is expected to be space-like. Furthermore, it turns out 
that most of the key results of physical interest j3(JII31) . such as the area increase law and general- 
ization of black hole mechanics, do not require the condition on the sign of CnQ{i)- It is therefore 
convenient to introduce a simpler and at the same time 'tighter' notion, that of a dynamical 
horizon, which is better suited to analyze how black holes grow in exact general relativity \'M\ : 

Definition 5: A smooth, three-dimensional, space-like sub-manifold (possibly with boundary) H 
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of space-time is said to be a dynamical horizon (DH) if it can be foliated by a family of closed 
2-manifolds such that 

1. on each leaf S the expansion O(^) of one null normal vanishes; and 

2. the expansion Q(n) of the other null normal n° is negative. 

Note first that, like FOTHs, dynamical horizons are 'space-time notions', defined quasi-locally. 
They are not defined relative to a space-like surface as was the case with Hawking's apparent 
horizons nor do they make any reference to infinity as is the case with event horizons. In partic- 
ular, they are well-defined also in the spatially compact context. Being quasi-local, they are not 
teleological. Next, let us spell out the relation between FOTHs and DHs. A space-like FOTH is a 
DH on which the additional condition £„ ^(e) < holds. Similarly, a DH satisfying £„ ^{i) < 
is a space-like FOTH. Thus, while neither definition implies the other, the two are closely related. 
The advantage of Definition [S] is that it refers only to the intrinsic structure of , without any 
conditions on the evolution of fields in directions transverse to H. Therefore, it is easier to verify 
in numerical simulations. More importantly, as we will see, this feature makes it natural to analyze 
the structure of H using only the constraint (or initial value) equations on it. This analysis will 
lead to a wealth of information on black hole dynamics. Reciprocally, Definition 0] has the advan- 
tage that, since it permits H to be space-like or null, it is better suited to analyze the transition 
to equilibrium |31| . 

A DH which is also a FOTH will be referred to as a space-like future outer horizon (SFOTH). 
To fully capture the physical notion of a dynamical black hole, one should require both sets of 
conditions, i.e., restrict oneself to SFOTHs. For, stationary black holes admit FOTHS and there 
exist space-times 156 which admit dynamical horizons but no trapped surfaces; neither can be 
regarded as containing a dynamical black hole. However, it is important to keep track of precisely 
which assumptions are needed to establish specific results. Most of the results reported in this 
review require only those conditions which are satisfied on DHs. This fact may well play a role 
in conceptual issues that arise while generalizing black hole thermodynamics to non-equilibrium 
situations^. 



2.2.2 Examples 

Let us begin with the simplest examples of space-times admitting DHs (and SFOTHs). These are 
provided by the spherically symmetric solution to Einstein's equations with a null fluid as source, 
the Vaidya metric |179lll37lll86) . (Further details and the inclusion of a cosmological constant are 
discussed in EBO J^^t as the Schwarzschild-Kruskal solution provides a great deal of intuition 
for general static black holes, the Vaidya metric furnishes some of the much needed intuition 
in the dynamical regime by bringing out the key differences between the static and dynamical 
situations. However, one should bear in mind that both Schwarzschild and Vaidya black holes are 
the simplest examples and certain aspects of geometry can be much more complicated in more 
general situations. The 4-metric of the Vaidya space-time is given by 

gab = - - ^^^7^) VawVfcW -t- 2V(„i;Vb)r + (Va6'V;,6' + sin^ eVa(p^b4>) , (U) 

where M{v) is any smooth, non-decreasing function of v. Thus, (w, r, 6, (/>) are the ingoing Eddington- 
Finkelstein coordinates. This is a solution of Einstein's equations, the stress-energy tensor Tab being 

■^Indeed, the situation is similar for black holes in equilibrium. While it is physically reasonable to restrict 
oneself to IHs, most results require only the WIH boundary conditions. The distinction can be important in certain 
applications, e.g., in finding boundary conditions on the quasi-equilibrium initial data at inner horizons. 
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given by 

Ta,^^VaVV,V, (12) 

where M = dM/dv. Clearly, Tab satisfies the dominant energy condition if M > 0, and vanishes if 
and only if M = 0. Of special interest to us are the cases illustrated in Figure^] AI{v) is non-zero 
until a certain finite retarded time, say w = 0, and then grows monotonically, either reaching an 
asymptotic value Mq as v tends to infinity (panel a), or, reaching this value at a finite retarded 
time, say v = vq, and then remaining constant (panel b). In either case, the space-time region 
u < is flat. 

Let us focus our attention on the metric 2-spheres, which are all given by w = const, and 
r = const.. It is easy to verify that the expansion of the outgoing null normal £° vanishes if and 
only if (v — const, and) r — 2GM{v). Thus, these are the only spherically symmetric marginally 
trapped surfaces MTSs. On each of them, the expansion 0(n) of the ingoing normal n° is negative. 
By inspection, the 3-metric on the world tube r — 2GM (v) of these MTSs has signature (-I-, -I-) 
when AI{v) is non-zero and (0,+,+) if M(v) is zero. Hence, in the left panel of Figure 0] the 
surface r = 2GM{v) is the DH H. In the right panel of Figure 21 the portion of this surface 
V < vo is the DH H, while the portion v > vq is a non-expanding horizon. (The general issue 
of transition of a DH to equilibrium is briefly discussed in Section |5|) Finally, note that at these 
MTSs, Cn9{^i) = — < 0. Hence in both cases, the DH is an SFOTH. Furthermore, in the case 
depicted in the right panel of Figure0]the entire surface r = 2GM{v) is a FOTH, part of which is 
dynamical and part null. 

This simple example also illustrates some interesting features which are absent in the stationary 
situations. First, by making explicit choices of M(v)^ one can plot the event horizon using, say, 
Mathematica |189| and show that they originate in the flat space-time region v < 0, in anticipation 
of the null fluid that is going to fall in after w = 0. The dynamical horizon, on the other hand, 
originates in the curved region of space-time, where the metric is time- dependent, and steadily 
expands until it reaches equilibrium. Finally, as Figures^jillustrate, the dynamical and event hori- 
zons can be well separated. Recall that in the equilibrium situation depicted by the Schwarzschild 
space-time, a spherically symmetric trapped surface passes through every point in the interior of 
the event horizon. In the dynamical situation depicted by the Vaidya space-time, they all lie in the 
interior of the DH. However, in both cases, the event horizon is the boundary of J~(X"*"). Thus, 
the numerous roles played by the event horizon in equilibrium situations get split in dynamical 
contexts, some taken up by the DH. 

What is the situation in a more general gravitational collapse? As indicated in the beginning of 
this section, the geometric structure can be much more subtle. Consider 3-manifolds t which are 
foliated by marginally trapped compact 2-surfaces S. We denote by the normal whose expansion 
vanishes. If the expansion of the other null normal n° is negative, r will be called a marginally 
trapped tube (MTT). If the tube t is space-like, it is a dynamical horizon. If it is time- like, it will be 
called time-like membrane. Since future directed causal curves can traverse time-like membranes 
in either direction, they are not good candidates to represent surfaces of black holes; therefore they 
are not referred to as horizons. 

In Vaidya metrics, there is precisely one MTT to which all three rotational Killing fields are 
tangential and this is the DH H. In the Oppenheimer-Volkoff dust collapse, however, the situation 
is just the opposite; the unique MTT on which each MTS S is spherical is time-like 180 47 . Thus 
we have a time-like membrane rather than a dynamical horizon. However, in this case the metric 
does not satisfy the smoothness conditions spelled out at the end of Section ^ and the global time- 
like character of r is an artifact of the lack of this smoothness. In the general perfect fluid spherical 
collapse, if the solution is smooth, one can show analytically that the spherical MTT is space-like 
at sufficiently late times, i.e., in a neighborhood of its intersection with the event horizon |1U2) . 
For the spherical scalar field collapse, numerical simulations show that, as in the Vaidya solutions, 
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Figure 4: Penrose diagrams of Schwarzschild-Vaidya metrics for which the mass function M{v) 
vanishes for v <Q The space-time metric is fiat in the past of v = (i.e., in the shaded 

region). In the left panel, as v tends to infinity, M vanishes and M tends to a constant value Mq. 
The space-like dynamical horizon H , the null event horizon E, and the time-like surface r — 2Mq 
(represented by the dashed line) all meet tangentially at In the right panel, for v > vq we 
have M = 0. Space-time in the future of v = Vq is isometric with a portion of the Schwarzschild 
space-time. The dynamical horizon H and the event horizon E meet tangentially at v — vq. In 
both figures, the event horizon originates in the shaded flat region, while the dynamical horizon 
exists only in the curved region. 
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the spherical MTT is space- hke everywhere |1()2| . Finally, the geometry of the numerically evolved 
MTTs has been examined in two types of non-spherical situations: the axi-symmetric collapse of a 
neutron star to a Kerr black hole and in the head-on collision of two non- rotating black holes |48) . 
In both cases, in the initial phase the MTT is neither space-like nor time-like all the way around its 
cross-sections S. However, it quickly becomes space-like and has a long space-like portion which 
approaches the event horizon. This portion is then a dynamical horizon. There are no hard results 
on what would happen in general, physically interesting situations. The current expectation is 
that the MTT of a numerically evolved black hole space-time which asymptotically approaches the 
event horizon will become space-like rather soon after its formation. Therefore most of the ongoing 
detailed work focuses on this portion, although basic analytical results are available also on how 
the time-like membranes evolve (see Appendix A of 

2.2.3 Uniqueness 

Even in the simplest, Vaidya example discussed above, our explicit calculations were restricted to 
spherically symmetric marginally trapped surfaces. Indeed, already in the case of the Schwarzschild 
space-time, very little is known analytically about non-spherically symmetric marginally trapped 
surfaces. It is then natural to ask if the Vaidya metric admits other, non-spherical dynamical 
horizons which also asymptote to the non-expanding one. Indeed, even if we restrict ourselves to 
the 3-manifold r = 2GM{v), can we find another foliation by non-spherical, marginally trapped 
surfaces which endows it with another dynamical horizon structure? These considerations illustrate 
that in general there are two uniqueness issues that must be addressed. 

First, in a general space-time {Ai,gab), can a space-like 3-manifold H be foliated by two distinct 
families of marginally trapped surfaces, each endowing it with the structure of a dynamical horizon? 
Using the maximum principle, one can show that this is not possible j93| . Thus, if H admits a 
dynamical horizon structure, it is unique. 

Second, we can ask the following question: How many DHs can a space-time admit? Since 
a space-time may contain several distinct black holes, there may well be several distinct DHs. 
The relevant question is if distinct DHs can exist within each connected component of the (space- 
time) trapped region. On this issue there are several technically different uniqueness results |27) . 
It is simplest to summarize them in terms of SFOTHs. First, if two non-intersecting SFOTHs 
H and H' become tangential to the same non-expanding horizon at a finite time (see the right 
panel in Figure , then they coincide (or one is contained in the other) . Physically, a more 
interesting possibility, associated with the late stages of collapse or mergers, is that H and H' 
become asymptotic to the event horizon. Again, they must coincide in this case. At present, one 
can not rule out the existence of more than one SFOTHs which asymptote to the event horizon if 
they intersect each other repeatedly. However, even if this were to occur, the two horizon geometries 
would be non-trivially constrained. In particular, none of the marginally trapped surfaces on H 
can lie entirely to the past of H' . 

A better control on uniqueness is perhaps the most important open issue in the basic framework 
for dynamical horizons and there is ongoing work to improve the existing results. Note however 
that all results of Sections |21 and El including the area increase law and the generalization of black 
hole mechanics, apply to all DHs (including the 'transient ones' which may not asymptote to the 
event horizon). This makes the framework much more useful in practice. 

The existing results also provide some new insights for numerical relativity |27| . First, suppose 
that a MTT t is generated by a foliation of a region of space-time by partial Cauchy surfaces Mt 
such that each MTS St is the outermost MTS in Mt. Then r can not be a time-like membrane. 
Note however that this does not imply that r is necessarily a dynamical horizon because r may 
be partially time-like and partially space-like on each of its marginally trapped surfaces S. The 
requirement that r be space-like - i.e., be a dynamical horizon - would restrict the choice of the 
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foliation Mt of space-time and reduce the unruly freedom in the choice of gauge conditions that 
numerical simulations currently face. A second result of interest to numerical relativity is the 
following. Let a space-time [M^gab] admit a DH H which asymptotes to the event horizon. Let 
Mo be any partial Cauchy surface in {M.,gab) which intersects H in one of the marginally trapped 
surfaces, say Sq. Then, So is the outermost marginally trapped surface - i.e., apparent horizon in 
the numerical relativity terminology - on Mq. 
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3 Area Increase Law 



As mentioned in the introduction, the dynamical horizon framework has led to a monotonicity 
formula governing the growth of black holes. In this section, we summarize this result. Our 
discussion is divided into three parts. The first spells out the strategy, the second presents a brief 
derivation of the basic formula, and the third is devoted to interpretational issues. 

3.1 Preliminaries 

The first law of black hole mechanics ^ tells us how the area of the black hole increases when it 
makes a transition from an initial equilibrium state to a nearby equilibrium state. The question 
we want to address is: Can one obtain an integral generalization to incorporate fully dynamical 
situations? Attractive as this possibility seems, one immediately encounters a serious conceptual 
and technical problem. For, the generalization requires, in particular, a precise notion of the flux of 
gravitational energy across the horizon. Already at null infinity, the expression of the gravitational 
energy flux is subtle: One needs the framework developed by Bondi, Sachs, Newman, Penrose, and 
others to introduce a viable, gauge invariant expression of this flux j54IIH;^l[l85| . In the strong field 
regime, there is no satisfactory generalization of this framework and, beyond perturbation theory, 
no viable, gauge invariant notion of the flux of gravitational energy across a general surface. 

Yet, there are at least two general considerations that suggest that something special may 
happen on DHs. Consider a stellar collapse leading to the formation of a black hole. At the 
end of the process, one has a black hole and, from general physical considerations, one expects 
that the energy in the final black hole should equal the total matter plus gravitational energy 
that fell across the horizon. Thus, at least the total integrated flux across the horizon should be 
well defined. Indeed, it should equal the depletion of the energy in the asymptotic region, i.e., 
the difference between the ADM energy and the energy radiated across future null infinity. The 
second consideration involves the Penrose inequality |157| introduced in Section Q Heuristically, 
the inequality leads us to think of the radius of a marginally trapped surface as a measure of the 
mass in its interior, whence one is led to conclude that the change in the area is due to influx 
of energy. Since a DH is foliated by marginally trapped surfaces, it is tempting to hope that 
something special may happen, enabling one to define the flux of energy and angular momentum 
across it. This hope is borne out. 

In the discussion of DHs fSectionsHland lT!^ we will use the following conventions (see FigureEl. 
The DH is denoted by H and marginally trapped surfaces that foliate it are referred to as cross- 
sections. The unit, time-like normal to H is denoted by r° with gabT'^'^T^ = ^1- The intrinsic 
metric and the extrinsic curvature of H are denoted by qab ■= gab + TaJb and Kab ■= Qa'^qb^^ cTd: 
respectively. D is the derivative operator on H compatible with gab, T^ab its Ricci tensor, and 7?. 
its scalar curvature. The unit space-like vector orthogonal to S and tangent to H is denoted by f". 
Quantities intrinsic to S are generally written with a tilde. Thus, the two-metric on S is qab and 
the extrinsic curvature oi S <Z H is, Kab '■= 5a '''DcTd] the derivative operator on {S, qab) is D and 
its Ricci tensor is TZab- Finally, we fix the rescaling freedom in the choice of null normals to cross- 
sections via £° := r°-|-f" and := r" — (so that £°na = —2). To keep the discussion reasonably 
focused, we will not consider gauge fields with non-zero charges on the horizon. Inclusion of these 
fields is not difficult but introduces a number of subtleties and complications which are irrelevant 
for numerical relativity and astrophysics. 
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Figure 5: H is a dynamical horizon, foliated by marginally trapped surfaces S. r" is the unit 
time-like normal to H and f" the unit space-like normal within H to the foliations. Although H is 
space-like, motions along f" can he regarded as 'time evolution with respect to observers at infinity '. 
In this respect, one can think of H as a hyperboloid in Minkowski space and S as the intersection 
of the hyperboloid with space-like planes. In the figure, H joins on to a weakly isolated horizon A 
with null normal 1°" at a cross-section So . 



3.2 Area increase law 

The qualitative result that the area as of cross-sections S increases monotonically on H follows 
immediately from the definition, 

K = q-^'D^rh = ^9"^V„(4 - n^) = -^e(„) > 0, (13) 

since 6(^) = and 6(„) < 0. Hence as increases monotonically in the direction of f". The 
non-trivial task is to obtain a quantitative formula for the amount of area increase. 

To obtain this formula, one simply uses the scalar and vector constraints satisfied by the Cauchy 
data {qab,Kab) on H: 

Hs:=n + K^- K'^'Kab = 16^GT,b t^t^, (14) 
:= Dk (K-' - Kq-') = SirGf"^ %q\, (15) 

where ^ 

fab = Tab - '^:^^9ab, (16) 

and Tab is the matter stress-energy tensor. The strategy is entirely straightforward: One fixes two 
cross-sections Si and S2 oi H, multiplies Hs and Hy with appropriate lapse and shift fields and 
integrates the result on a portion AiJ C H which is bounded by and S2. Somewhat surprisingly, 
if the cosmological constant is non-negative, the resulting area balance law also provides strong 
constraints on the topology of cross sections S. 

Specification of lapse N and shift N°- is equivalent to the specification of a vector field ^" = 
Nt°' -h N°- with respect to which energy-fiux across H is defined. The definition of a DH provides 
a preferred direction field, that along 1°-. Hence it is natural set = Nl°- = Nt"' -\- Nr". We 
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will begin with this choice and defer the possibihty of choosing more general vector fields until 
Section Ol 

The object of interest now is the flux of energy associated with = Ni"' across AiJ. We 
denote the flux of matter energy across AT? by ^matter- 

(?) — / ^ .1?acb .3t 



-^i?itter:= / Ta.^'^ed'V. (17) 

J AH 

By taking the appropriate combination of Equations H14f) and H15|l we obtain 

= TTTT^ / N{H + K''~ K'^'^Kab + - Xq-^^)) d'V. (18) 

167rG J^H 

Since H is foliated by compact 2-manifolds S, one can perform a 2 + 1 decomposition of various 
quantities on H . In particular, one first uses the Gauss-Codazzi equation to express TZ in terms 
of 7?,, Kab, and a total divergence. Then, one uses the identity 

r\Kab + Kab)^Qit)^Q (19) 

to simplify the expression. Finally one sets 

CJab = r^t" ^aib, C = 5"'^V,4. (20) 

(Note that aab is just the shear tensor since the expansion of vanishes.) Then, Equation p8() 
reduces to 

[ NTld^V ^WnG [ fabf^S!' (fV + ( N {\a\^ + 2\C\^) SV. (21) 

J AH J AH J AH 

To simplify this expression further, we now make a specific choice of the lapse N . We denote 
by R the area-radius function; thus R is constant on each S and satisfies 05 = At:R^. Since we 
already know that area increases monotonically, i? is a good coordinate on and using it the 
3-volume <PV on H can be decomposed as (PV = \dR\^^dR<P'V , where d denotes the gradient on 
H. Therefore calculations simplify if we choose 

N = \dR\ = Nr. (22) 

We will set NrI'^ = ^(^ny Then, the integral on the left side of Equation (|21(l becomes 

I NnTZd^V ^ dRind^V ^I{R2- Ri), (23) 

J AH JRi J 

where Ri and R2 are the (geometrical) radii of Si and S2, and 2 is the Gauss-Bonnet topological 
invariant of the cross-sections S. Substituting back in Equation H21|l one obtains 



I {R2 - Ri) = WnG [ (Tab - J^9ab) ?"6% d^V + f Nr {\<jf + 2\(f) d^V. 

J AH V ^''^^ J J AH 



(24) 



This is the general expression relating the change in area to fluxes across Aff . Let us consider its 
ramifications in the three cases, A being positive, zero, or negative: 

• If A > 0, the right side is positive definite whence the Gauss-Bonnet invariant X is positive 
definite, and the topology of the cross-sections 5" of the DH is necessarily that of S"^ . 



23 



• If A = 0, then 5 is either spherical or toroidal. The toroidal case is exceptional: If it 
occurs, the matter and the gravitational energy flux across H vanishes (see Section \'6.'i\ . 
the metric qab is flat, CnQ(i) = (so H can not be a FOTH), and CiQi^^) = 0. In view 
of these highly restrictive conditions, toroidal DHs appear to be unrelated to the toroidal 
topology of cross-sections of the event horizon discussed by Shapiro, Teukolsky, Winicour, 
and others |121l 11671 In the generic spherical case, the area balance law H24II becomes 

^{R. - R.) = d^V + ^ J^Nn + 2ICP) d^V. (25) 

• If A < 0, there is no control on the sign of the right hand side of Equation (|24|) . Hence, 
a priori any topology is permissible. Stationary solutions with quite general topologies are 
known for black holes which are asymptotically locally anti-de Sitter. Event horizons of these 
solutions are the potential asymptotic states of these DHs in the distant future. 

For simplicity, the remainder of our discussion of DHs will be focused on the zero cosmological 
constant case with 2-sphere topology. 

3.3 Energy flux due to gravitational waves 

Let us interpret the various terms appearing in the area balance law ()25f) . 

The left side of this equation provides us with the change in the horizon radius caused by the 
dynamical process under consideration. Since the expansion Q(f,) vanishes, this is also the change 
in the Hawking mass as one moves from the cross section 5*1 to 52- The first integral on the right 
side of this equation is the flux .?"^tter of matter energy associated with the vector field . The 

second term is purely geometrical and accompanies the term representing the matter energy flux. 

( 

Hence it is interpreted as the flux Jgrav of S^^j^y energy carried by the gravitational radiation: 

^ifal ■■=T7^ f Nn{\a\' + 2|Cn d'V. (26) 

A priori, it is surprising that there should exist a meaningful expression for the gravitational 
energy flux in the strong held regime where gravitational waves can no longer be envisaged as 
ripples on a flat space-time. Therefore, it is important to subject this interpretation to viability 
criteria analogous to the 'standard' tests one uses to demonstrate the viability of the Bondi flux 
formula at null infinity. It is known that it passes most of these tests. However, to our knowledge, 
the status is still partially open on one of these criteria. The situation can be summarized as 
follows: 

Gauge invariance 

Since one did not have to introduce any structure, such as coordinates or tetrads, which is 
auxiliary to the problem, the expression is obviously gauge invariant. This is to be contrasted 
with definitions involving pseudo-tensors or background fields. 

Positivity 

The energy flux H26|l is manifestly non-negative. In the case of the Bondi flux, positivity 
played a key role in the early development of the gravitational radiation theory. It was per- 
haps the most convincing evidence that gravitational waves are not coordinate artifacts but 
carry physical energy. It is quite surprising that a simple, manifestly non-negative expression 
can exist in the strong field regime of DHs. One can of course apply our general strategy 
to any space-like 3-surface H, foliated by 2-spheres. However, if fl is not a DH, the sign of 
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the geometric terms in the integral over AiJ can not be controlled, not even when H lies in 
the black hole region and is foliated by trapped (rather than marginally trapped) surfaces S. 
Thus, the positivity of J^grav is a rather subtle property, not shared by 3-surfaces which are 
foliated by non-trapped surfaces, nor those which are foliated by trapped surfaces; one needs 
a foliation precisely by marginally trapped surfaces. The property is delicately matched to 
the definition of DHs [HJ ■ 

Locality 

All fields used in Equation l|26(l are defined by the local geometrical structures on cross- 
sections of H. This is a non-trivial property, shared also by the Bondi-flux formula. However, 
it is not shared in other contexts. For example, the proof of the positive energy theorem by 
Witten |188) provides a positive definite energy density on Cauchy surfaces. But since it 
is obtained by solving an elliptic equation with appropriate boundary conditions at infinity, 
this energy density is a highly non-local function of geometry. Locality of J-g^Jv enables one 
to associate it with the energy of gravitational waves instantaneously falling across any cross 
section S. 

Vanishing in spherical symmetry 

The fourth criterion is that the flux should vanish in presence of spherical symmetry. Sup- 
pose H is spherically symmetric. Then one can show that each cross-section of S must be 
spherically symmetric. Now, since the only spherically symmetric vector field and trace-free, 
second rank tensor field on a 2-sphere are the zero fields, aab — and C° = 0. 

Balance law 

The Bondi-Sachs energy flux also has the important property that there is a locally defined 
notion of the Bondi energy E{C) associated with any 2-sphere cross-section C of future null 
infinity, and the difference E(Ci) — E{C2) equals the Bondi-Sachs flux through the portion 
of null infinity bounded by C2 and Ci. Does the expression (|26ll share this property? The 
answer is in the affirmative: As noted in the beginning of this section, the integrated flux is 
precisely the difference between the locally defined Hawking mass associated with the cross- 
section. In Section we will extend these considerations to include angular momentum. 

Taken together, the properties discussed above provide a strong support in favor of the inter- 
pretation of Equation H26|l as the ^(jj) -energy flux carried by gravitational waves into the portion 
Ai? of the DH. Nonetheless, it is important to continue to think of new criteria and make sure 
that Equation H2t)|l passes these tests. For instance, in physically reasonable, stationary, vacuum 
solutions to Einstein's equations, one would expect that the flux should vanish. However, on DHs 
the area must increase. Thus, one is led to conjecture that these space-times do not admit DHs. 
While special cases of this conjecture have been proved, a general proof is still lacking. Situation 
is similar for non-spherical DHs in spherically symmetric space-times. 

We will conclude this section with two remarks: 

• The presence of the shear term jap in the integrand of the flux formula (|26|l seems natural 
from one's expectations based on perturbation theory at the event horizon of the Kerr fam- 
ily |lU8ll65j . But the term |Cp is new and can arise only because H is space- like rather than 
null: On a null surface, the analogous term vanishes identically. To bring out this point, one 
can consider a more general case and allow the cross-sections S to lie on a horizon which is 
partially null and partially space-like. Then, using a 2 + 2 formulation |117j one can show 
that flux on the null portion is given entirely by the term |crp |28| . However, on the space- 
like portion, the term |Cp does not vanish in general. Indeed, on a DH, it cannot vanish in 
presence of rotation: The angular momentum is given by the integral of Ca^°', where (p" is 
the rotational symmetry. 
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• The flux refers to a specific vector field ^"^-j and measures the change in the Hawking mass 
associated with the cross-sections. However, this is not a good measure of the mass in 
presence of angular momentum (see, e.g., |34) for numerical simulations). Generalization of 
the balance law to include angular momentum is discussed in Section 14.21 
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4 Black Hole Mechanics 



As mentioned in the introduction, the discovery of the laws of black hole mechanics has led to 
fundamental insights in both classical and quantum gravity. In this section we discuss how the 
standard framework tied to stationary space-times can be extended using WIHs and DHs. 

4.1 Mechanics of weakly isolated horizons 

The isolated horizon framework has not only extended black hole mechanics, but it has also led 
to a deeper insight into the 'origin' of the laws of black hole mechanics. In this section we will 
summarize these developments using WIHs. Along the way we shall also obtain formulas for the 
mass and angular momentum of a WIH. For simplicity, in the main part of the discussion, we 
will restrict ourselves to type II (i.e., axi-symmetric) WIHs on which all matter fields vanish. 
Generalizations including various types of matter field can be found in |19l 1261 [T51 1751 1751 177| . 

4.1.1 The zeroth law 

The zeroth law of thermodynamics says that the temperature of a system in thermodynamic 
equilibrium is constant. Its counterpart for black hole mechanics says that surface gravity of a 
weakly isolated horizon is constant. This result is non-trivial because the horizon geometry is 
not assumed to be spherically symmetric; the result holds even when the horizon itself is highly 
distorted so long as it is in equilibrium. It is established as follows. 

Recall from Section l!j . 1 . 31 that the notion of surface gravity is tied to the choice of a null normal 
of the isolated horizon: := i'^uja- Now, using Equation in Definition |21 (of WIHs) , we 
obtain: 

CeL^a = 0. (27) 
Next, recall from Equation (O that the curl of uja is related to the imaginary part of ^'2: 

dw = 2(Im*2)e (28) 

where e is the natural area 2-form on A satisfying Cie = and £ ■ e = 0. Hence we conclude 
£ ■ dto — which in turn implies that is constant on the horizon: 

= CiLu = d{£ ■ uj) = dKi^ey (29) 

This completes the proof of the zeroth law. As the argument shows, given an NEH, the main 
condition jSJ in the definition of a WIH is equivalent to constancy of surface gravity. Note that 
no restriction has been imposed on ^2 which determines the mass and angular momentum mul- 
tipoles [22 : as emphasized above, the zeroth law holds even if the WIH is highly distorted and 
rapidly rotating. 

If electromagnetic fields are included, one can also show that the electric potential is constant 
on the horizon [221 ■ Finally, there is an interesting interplay between the zeroth law the action 
principle. Let us restrict ourselves to space-times which admit a non expanding horizon as inner 
boundary. Then the standard Palatini action principle is not well defined because the variation 
produces a non-vanishing surface term at the horizon. The necessary and sufficient condition for 
this surface term to vanish is precisely that the gravitational (and the electromagnetic) zeroth 
laws hold |26| . Consequently, the standard action principle is well-defined if inner boundaries are 
WIHs. 



27 



4.1.2 Phase space, symplectic structure, and angular momentum 

In field theories, conserved quantities such as energy and angular momentum can be universally 
defined via a Hamiltonian framework: they are the numerical values of Hamiltonians generating 
canonical transformations corresponding to time translation and rotation symmetries. In absence 
of inner boundaries, it is this procedure that first led to the notion of the ADM energy and angular 
momentum at spatial infinity [7]. At null infinity, it can also be used to define fiuxes of Bondi 
energy and angular momentum across regions of I"*" 'SS , and values of these quantities associated 
with any cross-section of J"*" [THlll85| . 

This procedure can be extended to allow inner boundaries which are WIHs. The first ingredient 
required for a Hamiltonian framework is, of course, a phase space. The appropriate phase space 
now consists of fields living in a region of space-time outside the black hole, satisfying suitable 
boundary conditions at infinity and horizon. Let Ai be the region of space-time that we are 
interested in. The boundary of M consists of four components: the time-like cylinder r at spatial 
infinity, two space-like surfaces Mi and M2 which are the future and past boundaries of M. , and 
an inner boundary A which is to be the WIH (see Figure inj. At infinity, all fields are assumed to 
satisfy the fall-off conditions needed to ensure asymptotic fiatness. To ensure that A is a type II 
horizon, one fixes a rotational vector field ip'^ on A and requires that physical fields on M. are such 
that the induced geometry on A is that of a type II horizon with (/?° as the rotational symmetry. 




Figure 6: The region of space-time M. under consideration has an internal boundary A and is 
bounded by two Cauchy surfaces Mi and M2 and the time-like cylinder Too at infinity. M is a 
Cauchy surface in M whose intersection with A is a spherical cross-section S and the intersection 
with Too is Soo, the sphere at infinity. 

Two Hamiltonian frameworks are available. The first uses a covariant phase space which 
consists of the solutions to field equations which satisfy the required boundary conditions 26 
Here the calculations are simplest if one uses a first order formalism for gravity, so that the 
basic gravitational variables are orthonormal tetrads and Lorentz connections. The second uses a 
canonical phase space consisting of initial data on a Cauchy slice M of M |55) . In the gravitational 
sector, this description is based on the standard ADM variables. Since the conceptual structure 
underlying the main calculation and the final results are the same, the details of the formalism are 
not important. For definiteness, in the main discussion, we will use the covariant phase space and 
indicate the technical modifications needed in the canonical picture at the end. 

The phase space T is naturally endowed with a (pre-)symplectic structure r2 - a closed 2-form 
(whose degenerate directions correspond to infinitesimal gauge motions). Given any two vector 
fields (i.e., infinitesimal variations) 61 and 62 on F, the action fl{6i, 62) of the symplectic 2-form on 
them provides a function on F. A vector field X on F is said to be a Hamiltonian vector field (i.e., 
to generate an infinitesimal canonical transformation) if and only if Cx^ — 0. Since the phase 
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space is topologically trivial, it follows that this condition holds if and only if there is a function 
H on T such that n{5,X) = 6H for all vector fields 5. The fimction H is called a Hamiltonian 
and X its Hamiltonian vector field; alternatively, H is said to generate the infinitesimal canonical 
transformation X. 

Since we are interested in energy and angular momentum, the infinitesimal canonical trans- 
formations X will correspond to time translations and rotations. As in any generally covariant 
theory, when the constraints are satisfied, values of Hamiltonians generating such diffeomorphisms 
can be expressed purely as surface terms. In the present case, the relevant surfaces are the sphere 
at infinity and the spherical section S — M C\ /S. of the horizon. Thus the numerical values of 
Hamiltonians now consist of two terms: a term at infinity and a term at the horizon. The terms 
at infinity reproduce the ADM formulas for energy and angular momentum. The terms at the 
horizon define the energy and angular momentum of the WIH. 

Let us begin with angular momentum (see jl5j for details). Consider a vector field 0" on 
M which satisfies the following boundary conditions: (i) At infinity, 0° coincides with a fixed 
rotational symmetry of the fiducial fiat metric; and, (ii) on A, it coincides with the vector field ip"" . 
Lie derivatives of physical fields along 0" define a vector field ATj^j on T. The question is whether 
this is an infinitesimal canonical transformation, i.e., a generator of the phase space symmetry. As 
indicated above, this is the case if and only if there exists a phase space function J^"^) satisfying: 

W^"^) (30) 

for all variations 5. If such a phase space function j''*' exists, it can be interpreted as the Hamil- 
tonian generating rotations. 

Now, a direct calculation [T^ shows that, in absence of gauge fields on A, one has: 

mx^^)) = 1 'e] -<5jit]M W^^'- (31) 

As expected, the expression for SJ^"^^ consists of two terms: a term at the horizon and a term 
at infinity. The term at infinity is the variation of the familiar ADM angular momentum ./adm 
associated with 0° . The surface integral at the horizon is interpreted as the variation of the horizon 
angular momentum J^. Since variations S are arbitrary, one can recover, up to additive constants, 
"^ADM ^^'^ from their variations, and these constants can be eliminated by requiring that both 
of these angular momenta should vanish in static axi-symmetric space-times. One then obtains: 

Ja := - ^ £(c^a¥'") = - ^ £/ I"i[*2] (32) 

where the function / on A is related to by 9a/ — eba^^- In the last step we have used 
Equation H28|l and performed an integration by parts. Equation 132|l is the expression of the 
horizon angular momentum. Note that all fields that enter this expression are local to the horizon 
and is not required to be a Killing field of the space-time metric even in a neighborhood of the 
horizon. Therefore, Ja can be calculated knowing only the horizon geometry of a type II horizon. 
We conclude our discussion of angular momentum with some comments: 

The Hamiltonian j(<^) 

It follows from Equation H31|l and (|32|) that the total Hamiltonian generating the rotation 
along (/)° is the difference between the ADM and the horizon angular momenta (apart from 
a sign which is an artifact of conventions). Thus, it can be interpreted as the angular 
momentum of physical fields in the space-time region A4 outside the black hole. 
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Relation to the Komar integral 

If 0° happens to be a space-time Killing field in a neighborhood of A, then Ja agrees with 
the Komar integral of 4'°' ^B]. If 0" is a global, space-time Killing field, then both J^dm 
as well as Ja agree with the Komar integral, whence the total Hamiltonian J^*^^ vanishes 
identically. Since the fields in the space-time region Ai are all axi-symmetric in this case, 
this is just what one would expect from the definition of j'-'f'l 

Ja for general Eixial fields cp"' 

If the vector field (j)"' is tangential to cross-sections of A, J^"^) continues to the generator of 
the canonical transformation corresponding to rotations along 0", even if its restriction to 
A does not agree with the axial symmetry ip"' of horizon geometries of our phase space fields. 
However, there is an infinity of such vector fields i^" and there is no physical reason to identify 
the surface term Ja arising from any one of them with the horizon angular momentum. 

Inclusion of gauge fields 

If non-trivial gauge fields are present at the horizon. Equation H32|l is incomplete. The horizon 
angular momentum Ja is still an integral over S; however it now contains an additional term 
involving the Maxwell field. Thus Ja contains not only the 'bare' angular momentum but 
also a contribution from its electromagnetic hair (see for details). 

Canonical phase space 

The conceptual part of the above discussion does not change if one uses the canonical phase 
space in place of the covariant. However, now the generator of the canonical transfor- 
mation corresponding to rotations has a volume term in addition to the two surface terms 
discussed above. However, on the constraint surface the volume term vanishes and the nu- 
merical value of the Hamiltonian reduces to the two surface terms discussed above. 



4.1.3 Energy, mass, and the first law 

To obtain an expression of the horizon energy, one has to find the Hamiltonian on T generating 
diffeomorphisms along a time translation symmetry i° on Ai. To qualify as a symmetry, at infinity 
i° must approach a fixed time translation of the fiducial flat metric. At the horizon, must be 
an infinitesimal symmetry of the type II horizon geometry. Thus, the restriction of to A should 
be a linear combination of a null normal and the axial symmetry vector ip"", 

= (33) 

where and the angular velocity fJ^t) are constants on A. 

However there is subtlety: Unlike in the angular momentum calculation where ((>"■ is required to 
approach a fixed rotational vector (p"- on A, the restriction of i° to A can not be a fixed vector field. 
For physical reasons, the constants and fl(^t) should be allowed to vary from one space-time 

to another; they are to be functions on phase space. For instance, physically one expects Q(^t) 
to vanish on the Schwarzschild horizon but not on a generic Kerr horizon. In the terminology of 
numerical relativity, unlike (/s", the time translation t° must be a live vector field. As we shall see 
shortly, this generality is essential also for mathematical reasons: without it, evolution along f° 
will not be Hamiltonian! 

At first sight, it may seem surprising that there exist choices of evolution vector fields t° for 
which no Hamiltonian exists. But in fact this phenomenon can also happen in the derivations of 
the ADM energy for asymptotically flat space-times in the absence of any black holes. Standard 
treatments usually consider only those t° that asymptote to the same unit time translation at 
inflnity for all space-times included in the phase space. However, if we drop this requirement and 
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choose a live t"" which approaches different asymptotic time-translations for different space-times, 
then in general there exists no Hamiltonian which generates diffeomorphisms along such a t". 
Thus, the requirement that the evolution be Hamiltonian restricts permissible t"' . This restriction 
can be traced back to the fact that there is a fixed fiducial flat metric at infinity. At the horizon, 
the situation is the opposite: The geometry is not fixed and this forces one to adapt to the 
space-time under consideration, i.e., to make it live. 

Apart from this important caveat, the calculation of the Hamiltonian is very similar to that for 
angular momentum. First, one evaluates the 1-form Y(t) on F whose action on any tangent vector 
field 5 is given by 

y(,)(<5):=f^(<5,X(0), (34) 

where X^^•^ is the vector field on T induced by diffeomorphisms along t°- . Once again, Yi^i){5) will 
consist of a surface term at infinity and a surface term at the horizon. A direct calculation yields 

y(t){5) = -^l-^aA - ^(t)6JA + 5E%^, (35) 

where K{t) '■— Bi^i t)£°'uja is the surface gravity associated with the restriction of i° to A, oa is 
the area of A, and i?^*DM ADM energy associated with The first two terms in the right 

hand side of this equation are associated with the horizon, while the £'adm term is associated 
with an integral at infinity. Since the term at infinity gives the ADM energy, it is natural to 
hope that terms at the horizon will give the horizon energy. However, at this point, we see an 
important difference from the angular momentum calculation. Recall that the right hand side of 
Equation (|31|) is an exact variation which means that J^"^"^ is well defined. However, the right hand 
side of Equation H35|l is not guaranteed to be an exact variation; in other words, A'(f) need not be 
a Hamiltonian vector field in phase space. It is Hamiltonian if and only if there is a phase space 
function i?^'' ~ the would be energy of the WIH - satisfying 

5E^l^ = ^SaA + n^t)SJA- (36) 

In particular, this condition implies that, of the infinite number of coordinates in phase space, £'^, 
and ri(t) can depend only on two: and Ja- 

Let us analyze Equation H36|l . Clearly, a necessary condition for existence of i?^-* is just the 
integrability requirement 

Since K(t) and Q.(t) are determined by Equation H37|l is a constraint on the restriction to the 
horizon of the time evolution vector field A vector field for which E"^^ exists is called a 
permissible time evolution vector field. Since Equation H36|l is precisely the first law of black hole 
mechanics, is permissible if and only if the first law holds. Thus the first law is the necessary 
and sufficient condition that the evolution generated by t"" is Hamiltonian! 

There are infinitely many permissible vector fields t°. To construct them, one can start with 
a suitably regular function kq of a a and Ja, find Bj^) so that Kjt) = kq, solve Equation 137|l to 
obtain ^{t), and find a permissible t° with i° — i^f" — ^(t)ip°' on A Each permissible i° 

defines a horizon energy e'^^ and provides a first law H36(l . A question naturally arises: Can one 
select a preferred tg or, alternatively, a canonical function Ko(aA, ^/a)? Now, thanks to the no-hair 
theorems, we know that for each choice of {a a, Ja), there is precisely one stationary black hole in 
vacuum general relativity: the Kerr solution. So, it is natural to set Ko{aA, Ja) = '^Kerr(aA, >/a), 
or, more explicitly, 

.o(«A,JA)^ ^^f^^^';^\^, , (38) 
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where Ra is the area radius of the horizon, Ra = (aA/47r)^/^. Via Equation H87|l . this choice then 
leads to 

^(t) = f^Kcrr(aA, J a) = 2GRa ~' ^^^^^ 

The associated horizon energy is then: 

e1"^ = ^^^/Ri+^G^- (40) 
This canonical horizon energy is called the horizon mass: 

Ma := e1"^ (41) 

Note that, its dependence on the horizon area and angular momentum is the same as that in the 
Kerr space-time. Although the final expression is so simple, it is important to keep in mind that this 
is not just a postulate. Rather, this result is derived using a systematic Hamiltonian framework, 
following the same overall procedure that leads to the definition of the ADM 4-momentum at 
spatial infinity. Finally, note that the quantities which enter the first law refer just to physical 
fields on the horizon; one does not have to go back and forth between the horizon and infinity. 
We will conclude with three remarks: 

Relation to the ADM and Bondi energy 

Under certain physically reasonable assumptions on the behavior of fields near future time- 
like infinity i+, one can argue that, if the WIH extends all the way to «+, then the difference 
-^^ADM — Ma equals the energy radiated across future null-infinity |14| . Thus, as one would 
expect physically, Ma is the mass that is left over after all the gravitational radiation has 
left the system. 

Horizon angular momentum and mass 

To obtain a well-defined action principle and Hamiltonian framework, it is essential to work 
with WIHs. However, the final expressions H^^2|l and (|4()|l of the horizon angular momentum 
and mass do not refer to the preferred null normals [£] used in the transition from an NEH 
to a WIH. Therefore, the expressions can be used on any NEH. This fact is useful in the 
analysis of transition to equilibrium fSection l4.3() and numerical relativity (Section l5.1|l . 

Generalizations of the first law 

The derivation of the first law given here can be extended to allow the presence of matter 
fields at the horizon |2H1 If gauge fields are present, the expression of the angular 

momentum has an extra term and the first law (|36|) also acquires the familiar extra term 
'$(5(5', representing work done on the horizon in increasing its charge. Again, all quantities 
are defined locally on the horizon. The situation is similar in lower and higher jl33| 
space-time dimensions. However, a key difference arises in the definition of the horizon mass. 
Since the uniqueness theorems for stationary black holes fail to extend beyond the Einstein- 
Maxwell theory in four space-time dimensions, it is no longer possible to assign a canonical 
mass to the horizon. However, as we will see in Section |H| the ambiguity in the notion of 
the horizon mass can in fact be exploited to obtain new insights into the properties of black 
holes and solitons in these more general theories. 

4.2 Mechanics of dynamical horizons 

The variations S in the first law (|36|) represent infinitesimal changes in equilibrium states of horizon 
geometries. In the derivation of Section l^lTl these variations relate nearby but distinct space-times 
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in each of which the horizon is in equihbrium. Therefore Equation is interpreted as the first 
law in a passive form. Physically, it is perhaps the active form of the first law that is of more direct 
interest where a physical process^ such as the one depicted in the right panel of Figure ^ causes 
a transition from one equilibrium state to a nearby one. Such a law can be established in the 
dynamical horizon framework. In fact, one can consider fully non-equilibrium situations, allowing 
physical processes in a given space-time in which there is a finite - rather than an infinitesimal - 
change in the state of the horizon. This leads to an integral version of the first law. 

Our summary of the mechanics of DHs is divided in to three parts. In the first, we begin with 
some preliminaries on angular momentum. In the second, we extend the area balance law (|25() by 
allowing more general lapse and shift functions, which leads to the integral version of the first law. 
In the third, we introduce the notion of horizon mass. 



4.2.1 Angular momentum balance 

As one might expect, the angular momentum balance law results from the momentum con- 
straint (|15() on the DH H. Fix any vector field ip°- on H which is tangential to all the cross-sections 
S of H, contract both sides of Equation H15|l with tf'^, and integrate the resulting equation over 
the region AH to obtain^ 

It is natural to identify the surface integrals with the generalized angular momentum Jg associated 
with cross-sections S and set 

where the overall sign ensures compatibility with conventions normally used in the asymptotically 
flat context, and where we have introduced an angular momentum density j"^ := —Kabf'^f^ for 
later convenience. The term 'generalized' emphasizes the fact that the vector field (/j" need not be 
an axial Killing field even on S; it only has to be tangential to our cross-sections. If ip"' happens 
to be the restriction of a space-time Killing field to S, then Jg agrees with the Komar integral. 
If the pair {qab, Kat) is spherically symmetric on S, as one would expect, the angular momenta 
associated with the rotational Killing fields vanish. 

Equation is a balance law; the right side provides expressions of fluxes of the generalized 
angular momentum across AH. The contributions due to matter and gravitational waves are 
cleanly separated and given by 

JLttcr^^ [ T,b?^^'d'V, J-^^^ = --!-/■ p-'C^q,,d'V, (44) 

with P'^'' = K'''' - Kq"^, so that 

As expected, if ip"' is a Killing vector of the three- metric qab, then the gravitational angular mo- 
mentum flux vanishes: jTg^av = 0- 

As with the area balance law, here we worked directly with the constraint equations rather 
than with a Hamiltonian framework. However, we could also have used, e.g., the standard ADM 

*Note that wc could replace fat with T^t because Qabf^V^ = 0. Thus the cosmological constant plays no role in 
this section. 
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phase space framework based on a Cauchy surface M with internal boundary S and the outer 
boundary at infinity. If 0° is a vector field on M which tends to Lp°- on S and to an asymptotic 
rotational symmetry at infinity, we can ask for the phase space function which generates the 
canonical transformation corresponding to the rotation generated by 0° . When the constraints are 
satisfied, as usual the value of this generating function is given by just surface terms. The term 
at infinity provides the total angular momentum and, as in Section 14.1.21 it is natural to interpret 
the surface term at S as the iy9-angular momentum of S. This term can be expressed in terms of 
the Cauchy data {qab,Kab) on M as 



where f" is the unit normal to S within M. However, since the right side involves the extrinsic 
curvature of M, in general the value of the integral is sensitive to the choice of M. Hence, the 
notion of the (^-angular momentum associated with an arbitrary cross-section is ambiguous. This 
ambiguity disappears if (/?° is divergence-free on S. In particular, in this case, one has J'^ — '^S' 
Thus, although the balance law H42|l holds for more general vector fields (/s", it is robust only when 
is divergence-free on S. (These considerations shed some light on the interpretation of the field 
C" in the area balance law (|25() . For, the form of the right side of Equation H43|l implies that the 
field vanishes identically on S if and only if Jg vanishes for every divergence-free ip"' on S. In 
particular then, if the horizon has non-zero angular momentum, the ^-contribution to the energy 
flux can not vanish.) 

Finally, for Jg to be interpreted as 'the' angular momentum, ip"' has to be a symmetry. An 
obvious possibility is that it be a Killing field of qab on S. A more general scenario is discussed in 
Section IHl 

4.2.2 Integral form of the first law 

To obtain the area balance law, in Section IX^ we restricted ourselves to vector fields = Nui"", 
i.e., to lapse functions Nn = \dR\ and shifts iV° = NbJ^- We were then led to a conservation 
law for the Hawking mass. In the spherically symmetric context, the Hawking mass can be taken 
to be the physical mass of the horizon. However, as the Kerr space-time already illustrates, in 
presence of rotation this interpretation is physically incorrect. Therefore, although Equation (|25|) 
continues to dictate the dynamics of the Hawking mass even in presence of rotation, a more general 
procedure is needed to obtain physically interesting conservation laws in this case. In the case of 
WIHs, the first law incorporating rotations required us to consider suitable linear combinations of 

and the rotational symmetry field i^" on the horizon. In the same spirit, on DHs, one has to 
consider more general vector fields than more general choices of lapses and shifts. 

As on WIHs, one first restricts oneself to situations in which the metric qab on H admits a 
Killing field ip"" so that Jg can be unambiguously interpreted as the angular momentum associated 
with each S. In the case of a WIH, t° was given by t° = c£°' -\- VLtp"" and the freedom was in 
the choice of constants c and Jl. On a DH, one must allow the corresponding coefficients to be 
'time-dependent'. The simplest generalization is to choose, in place of vector fields 



where VL is an arbitrary function of R, and the lapse Nr is given by Nr — \dr\ for any function r 
of R. Note that one is free to rescale Nr and Q. by functions of R so that on each cross-section 
('instant of time') one has the same rescaling freedom as on a WIH. One can consider even more 
general lapse-shift pairs to allow, e.g., for differential rotation (see |31|'). 




(46) 



t'' iv^r + nip" = N, 



V?" + (iV,P + r!(p'^), 



(47) 
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Using in place of one obtains the following generalization of the area balance equa- 
tion IS]: 



r2 - ri 1 



2G SttG 



'2 J Si Jui Js J 



f Tab^'^t'd^V+j^ f Nr{\a\^ + 2\C\^) d^V--^ f np-'C^qabd^V. (48) 

Note that there is one balance equation for every vector field of the form H47|l : as in Section ITTl 
we have an infinite number of relations, now ensured by the constraint part of Einstein's equations. 
The right side of Equation (|48|) can be naturally interpreted as the flux T^^^ of the 'energy' 
associated with the vector field t° across AH. Hence, we can rewrite the equation as 

^AH^'^^^ + ^i i ^fd?V- i nfd^V- f dnj^fd^V]. (49) 



'S2 



/ ^dnifd^V] 
Jui Js J 



If 5*1 and S2 are only infinitcsimally separated, this integral equation reduces to the differential 
condition 



SttG I 2R dR 



Sag + f2 6Jg 

s 



= —Sas + nSJl (50) 

Thus, the infinitesimal form of Equation (|48(l is a familiar first law, provided [{l/2R){dr /dR)]{S) is 
identified as an effective surface gravity on the cross-section S. This identification can be motivated 
as follows. First, on a spherically symmetric DH, it is natural to choose r = R. Then the surface 
gravity reduces to 1/{2R), just as one would hope from one's experience with the Schwarzschild 
metric and more generally with static but possibly distorted horizons (See Appendix A of 14 ). 
Under the change R t-^ ^iR), we have Rr = (dr/dR) Rji, which is the natural generalization of the 
transformation property = ck^ of surface gravity of WIHs under the change £ ^ ci. Finally, 
Kr can also be regarded as a 2-sphere average of a geometrically defined surface gravity associated 
with certain vector fields on El '^'31'; hence the adjective 'effective'. 

To summarize. Equation H48|l represents an integral generalization of the first law of mechanics 
of weakly isolated horizons to dynamical situations in which the horizon is permitted to make a 
transition from a given state to one far away, not just nearby. The left side represents the fiux 
•^AH of the energy associated with the vector field t° , analogous to the fiux of Bondi energy across 
a portion of null infinity. A natural question therefore arises: Can one integrate this flux to obtain 
an energy E]^ which depends only on fields defined locally on the cross-section, as is possible at 
null infinity? As discussed in the next section, the answer is in the affirmative and the procedure 
leads to a canonical notion of horizon mass. 



4.2.3 Horizon mass 

In general relativity, the notion of energy always refers to a vector field. On DHs, the vector field 
is Therefore, to obtain an unambiguous notion of horizon mass, we need to make a canonical 
choice of t"' = Nri"' + fit/j", i.e., of functions TV,- and fl on H. As we saw in Section [4.1.31 on 
WIHs of 4-dimensional Einstein-Maxwell theory, the pair (oa, Ja) suffices to pick a canonical 
time translation field on A. The associated horizon energy E*^ is then interpreted as the mass 
M^. This suggests that the pair (as,Jg) be similarly used to make canonical choices Nj.o and 
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fi" on S. Thanks to the black hole uniqueness theorems of the 4-dimensional Einstein-Maxwell 
theory, this strategy is again viable. 

Recall that the horizon surface gravity and the horizon angular velocity in a Kerr solution can 
be expressed as a function only of the horizon radius R and angular momentum J: 

Given a cross section 5 of H, the idea is to consider the unique Kerr solution in which the horizon 
area is given by as and angular momentum by Js, and assign to S effective surface gravity Rs and 
angular velocity fls through 

lis ■= KKcriiRs, Js)i := ^KcrARSi J^) (52) 

Repeating this procedure on every cross-section, one obtains functions R^{R) and Vf'{R) on iJ, 
since J'^ is a function of R alone. The definition of the effective surface gravity then determines 
a function r*^ of R and hence iV^o uniquely. Thus, using Equation (|51|) . one can select a canonical 
vector field and Equation (|49() then provides a canonical balance law: 



The key question is whether this equation is integrable, i.e., if 



s 



^'Ih = eX - eX (54) 

for some which depends locally on fields defined on S. The answer is in the affirmative. 
Furthermore, the expression of E^° is remarkably simple and is identified with the horizon mass: 

Thus, on any cross-section 5, Mg is just the mass of the Kerr space-time which has horizon area 
as and angular momentum J^: As far as the mass is concerned, one can regard the DH as an 
evolution through 'a sequence of Kerr horizons'. The non-triviality of the result lies in the fact 
that, although this definition of mass is so 'elementary', thanks to the balance law H48|l it obeys a 
Bondi-type fiux formula, 

Ms, - Ms, - / Tahiti <fV + -i- / Nro (lap -I- 2|Cn (i'V (56) 

for a specific vector field — N^oi"' + il'-'ip"-, where each term on the right has a well-defined 
physical meaning. Thus, DHs admit a locally-defined notion of mass and an associated, canonical 
conservation law (|56|l . The availability of a mass formula also provides a canonical integral version 
of the first law through Equation (|48|) : 

Ms, - Ms, = ^-^-^ + rcif '^'^ ~ / "^'^ ^ fy^° f^" '^'^] ■ ^^^^ 

^\'^52 >^ Si J O'^ J S I 

The infinitesimal version of this equation yields the familiar first law 5M — {R/8TrG)5a + QSJ. 

On weakly isolated as well as dynamical horizons, area a and angular momentum J arise as the 
fundamental quantities and mass is expressed in terms of them. The fact that the horizon mass is 
the same function of a and J in both dynamical and equilibrium situations is extremely convenient 
for applications to numerical relativity [SI]- Conceptually, this simplicity is a direct consequence 
of the first law and the non-triviality lies in the existence of a balance equation H48|l , which makes 
it possible to integrate the first law. 
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4.3 Passage of dynamical horizons to equilibrium 



In physical situations, such as a gravitational collapse or black hole mergers, one expects the 
dynamical horizon to approach equilibrium at late times and become isolated. Because of back 
scattering, generically the approach is only asymptotic. However, the back scattering is generally 
quite weak and in simulations, within numerical errors, equilibrium is reached at finite times. The 
passage to equilibrium can be studied in detail in Vaidya solutions discussed in Section 12.2.21 
Moreover, in spherically symmetric examples such as these solutions, exact equilibrium can be 
reached at a finite time (see right panel of FigureQJ. The question then arises: In these situations, 
do various notions introduced on dynamical horizons go over smoothly to those introduced on 
WIHs? This issue has been analyzed only in a preliminary fashion 18 Ij . In this section we will 
summarize the known results. 

First, if the dynamical horizon is a FOTH, as the flux of matter and shear across H tends 
to zero, H becomes null and furthermore a non-expanding horizon. By a suitable choice of null 
normals, it can be made weakly isolated. Conditions under which it would also become an isolated 
horizon arc not well-understood. Fortunately, however, the final expressions of angular momentum 
and horizon mass refer only to that structure which is already available on non-expanding horizons 
(although, as we saw in Section r4.1l the underlying Hamiltonian framework does require the horizon 
to be weakly isolated Therefore, it is meaningful to ask if the angular momentum and 

mass defined on the DHs match with those defined on the non-expanding horizons. In the case 
when the approach to equilibrium is only asymptotic, it is rather straightforward to show that the 
answer is in the affirmative. 

In the case when the transition occurs at a finite time, the situation is somewhat subtle. 
First, we now have to deal with both regimes and the structures available in the two regimes are 
entirely different. Second, since the intrinsic metric becomes degenerate in the transition from 
the dynamical to isolated regimes, limits are rather delicate. In particular, the null vector field 
i°- = t"" + on H diverges, while n° = t° — f" tends to zero at the boundary. A priori therefore, 
it is not at all clear that angular momentum and mass would join smoothly if the transition occurs 
at a finite time. However, a detailed analysis shows that the two sets of notions in fact agree. 

More precisely, one has the following results. Let Q = H U A he a. C*''+^ 3-manifold (with 
fc > 2), topologically S'^ x K. as in the second Penrose diagram of Figure 0] Let the space-time 
metric gab in a neighborhood of Q be C'' . The part iJ of Q is assumed to have the structure of a 
DH and the part A of a non-expanding horizon. Finally, the pull-back qab of gat to Q is assumed 
to admit an axial Killing field (p"-. Then we have: 

• The angular momentum J*^ and the mass M defined in the two regimes agree on the boundary 
between H and A. 

• The vector field Iq defined on H and used in the definition of mass matches with a preferred 
vector field <q used to define mass on A. 

This agreement provides an independent support in favor of the strategy used to introduce the 
notion of mass in the two regimes. 
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5 Applications in Numerical Relativity 



By their very nature, numerical simulations of space-times are invariably tied to choices of coordi- 
nates, gauge conditions, dynamical variables, etc. Therefore, it is non-trivial to extract from them 
gauge invariant physics, especially in the strong curvature regions. Traditionally, the analytical 
infrastructure available for this purpose has been based on properties of the Kerr solution and its 
perturbations. However, a priori it is not clear if this intuition is reliable in the fully dynamical, 
strong curvature regime. On the numerical side, a number of significant advances have occurred 
in this area over the past few years. In particular, efficient algorithms have been introduced to 
find apparent horizons (see, e.g., (B^, 168. 177 ), black hole excision techniques have been success- 
fully implemented [731 a-nd the stability of numerical codes has steadily improved [HJ]. To take 
full advantage of these ongoing improvements, one must correspondingly 'upgrade' the analytical 
infra-structure so that one can extract physics more reliably and with greater accuracy. 

These considerations provided stimulus for a significant body of research at the interface of 
numerical relativity and the dynamical and isolated horizon frameworks. In this section, we will 
review the most important of these developments. Section 15.11 summarizes calculations of mass 
and angular momentum of black holes. Section 15.21 discusses applications to problems involving 
initial data. Specifically, we discuss the issue of constructing the 'quasi-cquilibrium initial data' and 
the calculation of the gravitational binding energy for a binary black hole problem. Section 15.31 
describes how one can calculate the source multipole moments for black holes, and Section 15.41 
presents a 'practical' approach for extracting gauge invariant waveforms. Throughout this section 
we assume that vacuum equations hold near horizons. 

5.1 Numerical computation of black hole mass and angular momentum 

As we saw in Section 0] the mechanics of IHs and DHs provides expressions of angular momentum 
and mass of the horizon. These expressions involve geometric quantities defined intrinsically on 
the IH A and DH H. Numerical simulations, on the other hand, deal with the 3-metric qab and 
extrinsic curvature Kab on (partial) Cauchy surfaces M . Therefore, the first task is to recast the 
formulas in terms of this Cauchy data. 

Simulations provide us with a foliation of space-time Ai by partial Cauchy surfaces M, each 
of which has a marginally trapped 2-surface S as (a connected component of) its inner boundary. 
The world tube of these 2-surfaces is a candidate for a DH or an IH. If it is space-like, it is a DH 
H and if it is null (or, equivalently, if the shear aab of the outward null normal to S is zero) it is 
a WIH A. The situation is depicted in Figured It is rather simple to numerically verify if these 
restrictions are met. To calculate mass and angular momentum, one assumes that the intrinsic 
2-metric on the cross-sections S admits a rotational Killing field (/s" (see, however, Section |H1 for 
weakening of this assumption). A rather general and convenient method, based on the notion of 
Killing transport, has been introduced and numerically implemented to explicitly find this vector 
field ip" 113 . 

Let us first suppose that, in a neighborhood of the cross-section S of interest, the world tube of 
marginally trapped surfaces constitutes an IH. Then the task is to recast Equation H32(l in terms 
of the Cauchy data [qab, Kab) on M. This task is also straightforward jHS] and one arrives at^: 



where i?° is the unit radial normal to S in AI . This formula is particularly convenient numerically 
since it involves the integral of a single component of the extrinsic curvature. 

^This formula has a different sign from that given in 1851 due to a difference in the sign convention in the definition 
of the extrinsic curvature. 




(58) 
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Figure 7: The world tube of apparent horizons and a Cauchy surface M intersect in a 2-sphere S . 
is the unit time-like normal to M and i?" is the unit space-like normal to S within M . 

Now consider the dynamical regime, i.e., assume that, in a smaU neighborhood of 5, the 
world tube of marginally trapped surfaces is a DH H. The angular momentum formula H46(l 
on DHs involves the Cauchy data on H. However, it is easy to show 31 that it equals the 
expression (158(1 involving Cauchy data on M . Thus, Equation \5f^) is in fact applicable in both 
the isolated and dynamical regimes. The availability of a single formula is extremely convenient 
in numerical simulations. From now on, we will drop the subscript A and denote the angular 
momentum of S simply by j'^^^ : 



Jg is the intrinsic angular momentum (i.e., spin) of the black hole at the instant of 'time' 
represented by S. Note that Equation (|59l) differs from the standard formula for the ADM angular 
momentum Jadm only in that the integral is over the apparent horizon instead of the sphere at 
infinity. Finally, we emphasize that one only assumes that ^p"^ is a rotational Killing field of the 
intrinsic 2-metric on S\ it does not have to extend to a Killing field even to a neighborhood of 
S. In fact, the expressions (|32l) and (|46() on IHs and DHs are meaningful if (p" is replaced by any 
vector field 0" tangential to S\ on A the expression is conserved and on H it satisfies a balance law. 
Furthermore, they both equal Equation H59|l if (^^ is divergence-free; it has to Lie-drag only the area 
2-form (rather than the intrinsic metric) on S. However, since S admits an infinity of divergence- 
free vector fields (^°, there is no a priori reason to interpret j'g'^ as the 'angular momentum' of the 
black hole. For this interpretation to be meaningful, must be chosen 'canonically'. The most 
unambiguous way to achieve this is to require that it be a Killing field of the intrinsic geometry of 
5*. However, in Section |S1 we will introduce a candidate vector field which could be used in more 
general situations which are only 'near axi-symmetry'. 

Having calculated J5, it is easy to evaluate the mass Ms using Equation H4U|I : 



This method for calculating Jg and Mg has already been used in numerical simulations, particularly 
in the analysis of the gravitational collapse of a neutron star to a black hole 34, ,48^. It is found 
that this method is numerically more accurate than the other commonly used methods. 




(59) 




(60) 
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Comparison with other methods 

The approach summarized above has the key advantage that it is rooted in Haniihonian 
methods which can be universally used to define conserved quantities. In particular, it is a 
natural extension of the approach used to arrive at the ADM and Bondi-Sachs quantities 
at spatial and null infinity |71 1181 ITH^ . From more 'practical' considerations, it has three 
important features: 

• The procedure does not presuppose that the horizon geometry is precisely that of the 
Kerr horizon. 

• It is coordinate independent. 

• It only requires data that is intrinsic to the apparent horizon. 

None of the commonly used alternatives share all three of these features. 

Before the availability of the IH and DH frameworks, standard procedures of calculating 
angular momentum were based on properties of the Kerr geometry. The motivation comes 
from the common belief, based on black-hole uniqueness theorems, that a black hole created 
in a violent event radiates away all its higher multipole moments and as it settles down, its 
near horizon geometry can be approximated by that of the Kerr solution. The strategy then 
is to identify the geometry of S with that of a suitable member of the Kerr family and read 
off the corresponding angular momentum and mass parameters. 

The very considerations that lead one to this strategy also show that it is not suitable in 
the dynamical regime where the horizon may be distorted and not well-approximated by 
any Kerr horizon. For horizons which have very nearly reached equilibrium, the strategy is 
physically well motivated. However, even in this case, one has to find a way to match the 
horizon of the numerical simulation with that of a specific member in the Kerr family. This 
is non-trivial because the coordinate system used in the given simulation will, generically, 
not bear any relation to any of the standard coordinate systems used to describe the Kerr 
solution. Thus, one cannot just look at, say, a metric component to extract mass and angular 
momentum. 

A semi-heuristic but most-commonly used procedure is the great circle method. It is based on 
an observation of the properties of the Kerr horizon made by Smarr jlfi9| using Kerr-Schild 
coordinates. Let Lg be the length of the equator and Lp the length of a polar meridian 
on the Kerr horizon, where the equator is the coordinate great circle of maximum proper 
length and a polar meridian is a great circle of minimum proper length. Define a distortion 
parameter 5 as 6 = (L^ — Lp)/ L^. The knowledge of 5, together with one other quantity such 
as the area, Lo, or Lp, is sufficient to find the parameters m and a of the Kerr geometry. 
However, difficulties arise when one wishes to use these ideas to calculate M and a for a 
general apparent horizon S. For, notions such as great circles, equator or polar meridian are 
all highly coordinate dependent. Indeed, if we represent the standard two-metric on the Kerr 
horizon in different coordinates, the great circles in one coordinate system will not agree with 
great circles in the other system. Therefore, already for the Kerr horizon, two coordinate 
systems will lead to different answers for A'l and a! In certain specific situations where 
one has a good intuition about the coordinate system being used and the physical situation 
being modelled, this method can be useful as a quick way of estimating angular momentum. 
However, it has the conceptual drawback that it is not derived from a well-founded, general 
principle and the practical drawback that it suffers from too many ambiguities. Therefore it 
is inadequate as a general method. 

Problems associated with coordinate dependence can be satisfactorily resolved on axi-sym- 
metric horizons, even when the coordinate system used in the numerical code is not adapted 
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to the axial symmetry. The idea is to use the orbits of the KiUing vector as analogs of the 
lines of latitude on a metric two-sphere. The analog of the equator is then the orbit of the 
Killing vector which has maximum proper length. This defines Le in an invariant way. The 
north and south poles are the points where the Killing vector vanishes, and the analog of Lp 
is the length of a geodesic joining these two points. (Because of axial symmetry, all geodesies 
joining the poles will have the same length). This geodesic is necessarily perpendicular to 
the Killing vector. Hence one just needs to find the length of a curve joining the north and 
south poles which is everywhere perpendicular to the Killing orbits. With and Lp defined 
in this coordinate invariant way, one can follow the same procedure as in the great circle 
method to calculate the mass and angular momentum. This procedure has been named |85| 
the generalized great circle method. 

How does the generalized great circle method compare to that based on IHs and DHs? Since 
one must find the Killing vector, the first step is the same in the two cases. In the IH and DH 
method, one is then left simply with an integration of a component of the extrinsic curvature 
on the horizon. In the generalized great circle method, by contrast, one has to determine 
the orbit of the Killing vector with maximum length and also to calculate the length of a 
curve joining the poles which is everywhere orthogonal to the Killing orbits. Numerically, 
this requires more work and the numerical errors are at least as large as those in the IH- 
DH method. Thus, even if one ignores conceptual considerations involving the fundamental 
meaning of conserved quantities, and furthermore restricts oneself to the non-dynamical 
regime, the practical simplicity of the great circle method is lost when it is made coordinate 
invariant. To summarize, conceptually. Equations H59|l and (I60|) provide the fundamental 
definitions of angular momentum and mass, while the great circle method provides a quick 
way of estimating these quantities in suitable situations. By comparing with Equations (|59l) 
and (|6()(l , one can calculate errors and sharpen intuition on the reliability of the great circle 
method. 

A completely different approach to finding the mass and angular momentum of a black hole 
in a numerical solution is to use the concept of a Killing horizon. Assume the existence of 
Killing vectors in the neighborhood of the horizon so that mass and angular momentum are 
defined as the appropriate Komar integrals. This method is coordinate independent and does 
not assume, at least for angular momentum, that the near horizon geometry is isometric with 
the Kerr geometry. But it has two disadvantages. First, since the Komar integral can involve 
derivatives of the Killing field away from the horizon, one has to find the Killing fields in a 
neighborhood of the horizon. Second, existence of such a stationary Killing vector is a strong 
assumption; it will not be satisfied in the dynamical regime. Even when the Killing field 
exists, computationally it is much more expensive to find it in a neighborhood of the horizon 
rather than the horizon itself. Finally, it is not a priori clear how the stationary Killing vector 
is to be normalized if it is only known in a neighborhood of the horizon. In a precise sense, the 
isolated horizon framework extracts just the minimum amount of information from a Killing 
horizon in order to carry out the Hamiltonian analysis and define conserved quantities by 
by-passing these obstacles. 

5.2 Initial data 

In order to accurately calculate the gravitational waveforms for, say, the coalescence of binary 
black holes, one must begin with astrophysically relevant initial data. While there has been some 
progress on the construction of such data, the general problem is yet to be solved (see, e.g., [Slj for 
a recent review). However, there has been notable progress in the case of black holes which are in 
quasi-equilibrium. Physically, this situation arises when the black holes are sufficiently far apart for 
their orbits to be quasi-periodic. In his section, we will first summarize this development [T^IT^ 
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and then review an approach to the calculation of the binding energy in the initial data. The 
discussion is based on the IH framework. 

5.2.1 Boundary conditions at the inner boundary 

Consider the problem of constructing initial data {qab, Kab) on M , representing a binary black hole 
system. This problem has two distinct aspects. The first has to do with the inner boundaries. A 
natural avenue to handle the singularities is to 'excise' the region contained in each black hole (see, 
e.g., pil73p. However, this procedure creates inner boundaries on M and one must specify suitable 
boundary conditions there. The boundary conditions should be appropriate for the elliptic system 
of constraint equations, and they should also capture the idea that the excised region represents 
black holes which have certain physical parameters and which are in quasi-equilibrium at the 
'instant' represented by M . The second aspect of the problem is the choice of the free data in 
the bulk. To be of physical interest, not only must the free data satisfy the appropriate boundary 
conditions, but the values in the bulk must also have certain properties. For example, we might 
want the black holes to move in approximately circular orbits around each other, and require that 
there be only a minimal amount of spurious gravitational radiation in the initial data. The black 
holes would then remain approximately isolated for a sufficiently long time, and orbit around each 
other before finally coalescing. 

While a fully satisfactory method of prescribing such initial data is still lacking, there has 
been significant progress in recent years. When the black holes are far apart and moving on 
approximately circular orbits, one might expect the trajectory of the black holes to lie along an ap- 
proximate helical Killing vector [23 1^ llOiil 11041 0] . Using concepts from the helical Killing vector 
approximation and working in the conformal thin-sandwich decomposition of the initial data |192| . 
Cook has introduced the 'quasi-equilibrium' boundary conditions which require that each of the 
black holes be in instantaneous equilibrium [75|; see also |190l I74| for a similar approach. The 
relation between these quasi-equilibrium boundary conditions and the isolated horizon formalism 
has also been recently studied |124| . 

In this section, we consider only the first aspect mentioned above, namely the inner bound- 
aries. Physically, the quasi-equilibrium approximation ought to be valid for time intervals much 
smaller than other dynamical time scales in the problem, and the framework assumes only that 
the approximation holds infinitesimally 'off Af. So, in this section, the type II NEH A will be an 
infinitesimal world tube of apparent horizons. We assume that there is an axial symmetry vector 
(/?° on the horizon, although, as discussed in Section this assumption can be weakened. 

Depending on the degree of isolation one wants to impose on the individual black holes, the 
inner boundary may be taken to be the cross-section of either an NEH, WIH, or an IH. The strategy 
is to first start with an NEH and impose successively stronger conditions if necessary. Using the 
local geometry of intersections S' of M with A, one can easily calculate the area radius i?A, the 
angular momentum Ja given by Equation H59|l . the canonical surface gravity k = KKcrr(-RA, J a) 
given by Equation (|38|l . and angular velocity Vl — r2Korr(-RAj ^/a) given by Equation (jSHl- (Note 
that while a WIH structure is used to arrive at the expressions for Ja, and fi, the expressions 
themselves are unambiguously defined also on a NEH.) These considerations translate directly into 
restrictions on the shift vector /?" at the horizon. If, as in Section 14.1.31 one requires that the 
restriction of the evolution vector field i° to A be of the form 

e = Btl - nKcr-r-^" = BT'' + [BR'' - riKorr'/'"), (61) 

where — T" + i?° , then the values of the lapse and shift fields at the horizon are given by 

a = B, 13°- = aR" - f^KcrrVs"- (62) 
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We are free to choose the lapse a arbitrarily on S, but this fixes the shift The tangential part 
of /3° can, if desired, be chosen differently depending on the asymptotic value of (3°" at infinity 
which determines the angular velocity of inertial observers at infinity. 

Next, one imposes the condition that the infinitesimal world-tube of apparent horizons is an 
'instantaneous' non-expanding horizon. This requirement is equivalent to 

Ctqab = 0, C^Qab = (63) 

on A, which imply that the shear and expansion of vanish identically. These can be easily 
translated to conditions that the Cauchy data {qab,Kab) must satisfy at the horizon 1124) . 
These boundary conditions are equivalent to the quasi-equilibrium boundary conditions developed 
by Cook However, the isolated horizon formalism allows greater control on the physical 

parameters of the black hole. In particular, we are not restricted just to the co-rotational or 
irrotational cases. Furthermore, the use of IHs streamlines the procedure: Rather than adding 
conditions one by one as needed, one begins with a physically well-motivated notion of horizons in 
equilibrium and systematically derives its consequences. 

Up to this point, the considerations are general in the sense that they are not tied to a particular 
method of solving the initial value problem. However, for the quasi-equilibrium problem, it is the 
conformal thin-sandwich method |192II7()) that appears to be best suited. This approach is based on 
the conformal method [14411191] where we write the 3-metric as qab = V'^'Zafc- The free data consists 
of the conformal 3-metric Qabi its time derivative CtQab, the lapse a, and the trace of the extrinsic 
curvature K. Given this free data, the remaining quantities, namely the conformal factor and the 
shift, are determined by elliptic equations provided appropriate boundary conditions are specified 
for them on the horizon^. It turns out |124) that the horizon conditions 1)63)1 are well-tailored 
for this purpose. While the issue of existence and uniqueness of solutions using these boundary 
conditions has not been proven, it is often the case that numerical calculations are convergent and 
the resulting solutions are well behaved. Thus, these conditions might therefore be sufficient from 
a practical point of view. 

In the above discussion, the free data consisted of {qab,J~-tqab,C(,K) and one solved elliptic 
equations for {ip^P""). However, it is common to consider an enlarged initial value problem by 
taking LtK as part of the free data (usually set to zero) and solving an elliptic equation for a. We 
now need to prescribe an additional boundary condition for a. It turns out that this can be done 
by using WIHs, i.e., by bringing in surface gravity, which did not play any role so far. From the 
definition of surface gravity in Equations (0 IH} , it is clear that the expression for k will involve a 
time derivative; in particular, it turns out to involve the time derivative of a. It can be shown that 
by choosing Cta on S (e.g., by taking Cta = 0) and requiring surface gravity to be constant on S 
and equal to KKerr, one obtains a suitable boundary condition for a. (The freedom to choose freely 
the function Cta mirrors the fact that fixing surface gravity does not uniquely fix the rescaling 
freedom of the null normal.) Note that k is required to be constant only on S, not on A. To 
ask it to be constant on A would require CtK — 0, which in turn would restrict the second time 
derivative; this necessarily involves the evolution equations, and they are not part of the initial 
data scheme. 

One may imagine using the yet stronger notion of an IH, to completely fix the value of the 
lapse at the horizon. But this requires solving an elliptic equation on the horizon and the relevant 
elliptic operator has a large kernel [161 172) . Nonetheless, the class of initial data on which its 
inverse exists is infinite dimensional so that the method may be useful in practice. However, this 
condition would genuinely restrict the permissible initial data sets. In this sense, while the degree 

^Cook's boundary condition on the conformal factor tp (Equation (82) in 1721 ') is equivalent to @[e) =0 which 
(in the co-rotating case, or more generally, when the 2-metric on S is axi-symmetric) reduces to Ct4' = on S. 
The Yo et al. boundary condition on ip (Equation (48) of [TgDl ) is equivalent to Cj;4' = on S, where, however, the 
evolution vector field is obtained by a superposition of two Kerr-Schild data. 
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of isolation implied by the IH boundary condition is likely to be met in the asymptotic future, for 
quasi-equilibrium initial data it is too strong in general. It is the WIH boundary conditions that 
appear to be well-tailored for this application. 

Finally, using methods introduced by Dain |80| . a variation of the above procedure was recently 
introduced to establish the existence and uniqueness of solutions and to ensure that the conformal 
factor ip is everywhere positive 81 . One again imposes Equation H63|l . However, in place of Dirich- 
let boundary conditions (|62(l on the shift, one now imposes Neumann- type conditions on certain 
components of /3". This method is expected to be applicable all initial data constructions relying 
on the conformal method. Furthermore, the result might also be of practical use in numerical 
constructions to ensure that the codes converge to a well behaved solution. 

5.2.2 Binding energy 

For initial data representing a binary black hole system, the quantity Eh = Madm — Mi — M2 is 
called the effective binding energy, where Madm is the ADM mass, and Mi^2 are the individual 
masses of the two black holes. Heuristically, even in vacuum general relativity, one would expect 
Eh to have several components. First there is the analog of the Newtonian potential energy and 
the spin-spin interaction, both of which may be interpreted as contributing to the binding energy. 
But Eh also contains contributions from kinetic energy due to momentum and orbital angular 
momentum of black holes, and energy in the gravitational radiation in the initial data. It is only 
when these are negligible that Eh is a good measure of the physical binding energy. 

The first calculation of binding energy was made by Brill and Lindquist in such a context. They 
considered two non-spinning black holes initially at rest For large separations, they found 

that, in a certain mathematical sense, the leading contribution to binding energy comes just from 
the usual Newtonian gravitational potential. More recently, Dain [73j has extended this calculation 
to the case of black holes with spin and has shown that the spin-spin interaction energy is correctly 
incorporated in the same sense. 

In numerical relativity, the notion of binding energy has been used to locate sequences of 
quasi-circular orbits. The underlying heuristic idea is to minimize Eh with respect to the proper 
separation between the holes, keeping the physical parameters of the black holes fixed. The value 
of the separation which minimizes Eh provides an estimate of sizes of stable 'circular' orbits |71l 
I39II159| . One finds that these orbits do not exist if the orbital angular momentum is smaller than 
a critical value (which depends on other parameters) and uses this fact to approximately locate the 
'inner- most stable circular orbit' (ISCO). In another application, the binding energy has been used 
to compare different initial data sets which are meant to describe the same physical system. If the 
initial data sets have the same values of the black hole masses, angular momenta, linear momenta, 
orbital angular momenta, and relative separation, then any differences in Eh should be due only 
to the different radiation content. Therefore, minimization of Eh corresponds to minimization of 
the amount of radiation in the initial data jlSB) . 

In all these applications, it is important that the physical parameters of the black holes are cal- 
culated accurately. To illustrate the potential problems, let us return to the original Brill-Lindquist 
calculation jHSl- The topology of the spatial slice is R'^ with two points (called 'punctures') re- 
moved. These punctures do not represent curvature singularities. Rather, each of them represents 
a spatial infinity of an asymptotically fiat region which is hidden behind an apparent horizon. This 
is a generalization of the familiar Einstein-Rosen bridge in the maximally extended Schwarzschild 
solution. The black hole masses Mi and M2 are taken to be the ADM masses of the corresponding 
hidden asymptotic regions. (Similarly, in j79|. the angular momentum of each hole is defined to be 
the ADM angular momentum at the corresponding puncture.) Comparison between Eh and the 
Newtonian binding energy requires us to define the distance between the holes. This is taken to be 
the distance between the punctures in a fiducial flat background metric; the physical distance be- 
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tween the two punctures is infinite since they represent asymptotic ends of the spatial 3-manifold. 
Therefore, the sense in which one recovers the Newtonian binding energy as the leading term is 
physically rather obscure. 

Let us re-examine the procedure with an eye to extending it to a more general context. Let 
us begin with the definition of masses of individual holes which are taken to be the ADM masses 
in the respective asymptotic regions. How do we know that these are the physically appropriate 
quantities for calculating the potential energy? Furthermore, there exist initial data sets (e.g., 
Misner data |15Ulll5l] ) in which each black hole does not have separate asymptotic regions; there 
are only two common asymptotic regions connected by two wormholes. For these cases, the use of 
ADM quantities is clearly inadequate. The same limitations apply to the assignment of angular 
momentum. 

A natural way to resolve these conceptual issues is to let the horizons, rather than the punctures, 
represent black holes. Thus, in the spirit of the IH and DH frameworks, it is more appropriate to 
calculate the mass and angular momentum using expressions H60I I59|l which involve the geometry 
of the two apparent horizons. (This requires the apparent horizons to be axi-symmetric, but 
this limitation could be overcome following the procedure suggested in Section |H1) Similarly, the 
physical distance between the black holes should be the smallest proper distance between the two 
apparent horizons. To test the viability of this approach, one can repeat the original Brill-Lindquist 
calculation in the limit when the black holes are far apart |136) . One first approximately locates 
the apparent horizon, finds the proper distance d between them, and then calculates the horizon 
masses (and thereby Sb) as a power series in 1/d. The leading-order term does turn out to be 
the usual Newtonian gravitational potential energy but the higher order terms are now different 
from |63| . Similarly, it would be interesting to repeat this for the case of spinning black holes and 
recover the leading order term of |79) within this more physical paradigm using, say, the Bowen- 
York initial data. This result would re-enforce the physical ideas and the approach can then be 
used as a well defined method for calculating binding energy in more general situations. 



In the examples mentioned above, we focussed on mass and angular momentum; linear 
momentum was not considered. A meaningful definition for the linear momentum of a black 
hole would be useful for several problems: It would enable one to define the orbital angular 
momentum, locate quasi-circular orbits and enable a comparison between different initial 
data sets. However, just as angular momentum is associated with rotational symmetry, 
linear momentum ought to be associated with a space translational symmetry. In a general 
curved space-time, one does not have even an approximate notion of translations except at 
infinity. Thus, whenever one speaks of linear momentum of initial data (e.g., in the Bowen- 
York construction jHZ]) the quantity is taken to be the appropriate component of the ADM 
momentum at infinity. For reasons discussed above, it is more appropriate to associate black 
hole parameters with the respective horizons. So, a question naturally arises: Can one define 
linear momentum as an integral over the apparent horizon? If z°' is a vector field on M 
which, for some reason, can be regarded as an approximate translational symmetry, and S 
the apparent horizon, then one might naively define the linear momentum associated with 



Unfortunately, this definition is not satisfactory. For example, in the boosted Kerr-Schild 
metric, if one lets z"' to be the translational vector field of the background Minkowski space. 
Equation (|64|) does not yield the physically expected linear momentum of the boosted black 
hole. See 136 for further discussion and for other possible ways to define linear momentum. 
Further work is needed to address this question satisfactorily. 



Remark 




(64) 
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5.3 Black hole multipole moments 

Let us begin by considering the notion of source and field multipole moments in Newtonian gravity 
and in flat space electrodynamics. Field multipoles appear in the asymptotic expansions of the 
fields at infinity while the source multipoles are defined in terms of the mass or charge distribution of 
the source. These two sets of multipole moments are related to each other via field equations. The 
same is true in linearized general relativity |155l . Also, the well known quadrupole formula relates 
the rate of change of the quadrupole moment to the energy flux at infinity due to gravitational 
waves. 

The situation in exact, non-linear general relativity is not so simple. Using the geometric 
structure of the gravitational field near spatial infinity, the field multipoles for stationary space- 
times were studied by Geroch, Hansen, Beig, Simon, and others piTI 11071 1^ HTl BHl They 
found that, just as in electrodynamics, the gravitational field has two sets of multipoles: The 
mass multipoles M„ and the angular momentum multipoles J„. The knowledge of these multipole 
moments suffices to determine the space-time geometry in a neighborhood of spatial infinity |41l 
1401 142| . Thus, at least in the context of stationary space-times, the field multipole moments are 
well understood. However, in problems involving equations of motion, it is the source multipoles 
that are of more direct interest, ft is natural to ask if these can be defined for black holes. 

The answer is affirmative for black holes in equilibrium, which can be represented by isolated 
horizons. For simplicity, we will consider only type II (i.e., axisymmetric) , non-extremal isolated 
horizons in vacuum. The source multipoles are two sets M„ and J„ of numbers which provide a 
diffeomorphism invariant characterization of the horizon geometry. 

As before, let 5 be a cross-section of A. We denote the intrinsic Riemannian metric on it by 
(jab, the corresponding area 2-form by iab, and the derivative operator by Da- Since the horizon 
is of type II, there exists a vector field on S such that Cipijab = 0. The two points where (p"- 
vanishes are called the poles of S. The integral curves of t^" are natural candidates for the 'lines 
of latitude' on S, and the lines of longitude are the curves which connect the two poles and are 
orthogonal to </?°. This leads to an invariantly defined coordinate C £ [—1, 1] ~ the analog of the 
function cos 9 in usual spherical coordinates - defined by 

DaC = -^ha^\ (65) 

where R is the area radius of S. The freedom of adding a constant to ^ is removed by requiring 
<fg Qi — 0. With (f) G [0, 2tt) being an affine parameter along the 2-metric qab takes the canonical 
form 

qab - RHf-^DaCDbC + fDa4>Db^), (66) 

where / — ipa'p"' / ■ The only remaining freedom in the choice of coordinates is a rigid shift in (j). 
Thus, in any axi-symmetric geometry, there is an invariantly defined coordinate C, and multipole 
moments are defined using the Legendre polynomials in C, as weight functions. 

Recall from Section 12.1.31 that the invariant content in the geometry of an isolated horizon is 
coded in (the value of its area and) ^2- The real part of is proportional to the scalar curvature 
R of qab and captures distortions |98II88| . while the imaginary part of 5*2 yields the curl of uja and 
captures the angular momentum information. (The free function / in Equation (|66l) determines 
and is completely determined by the scalar curvature R ) Multipoles are constructed directly from 
\E'2. The angular momentum multipoles are defined as 

where Y^{C,) are the spherical harmonics. (Axi-symmetry ensures that y,™ vanish if m ^ 0.) 
The normalization factors have been chosen to ensure that the dimensions of multipoles are the 
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physically expected ones and Ji is the angular momentum of the isolated horizon. The mass 
multipoles are defined similarly as the moments of R: 

where the normalization factor has been chosen in order to ensure the correct physical dimensions 
and Mo = A/a- 

These multipoles have a number of physically desired properties: 

• The angular momentum monopole vanishes, Jo = 0. (This is because we are considering only 
smooth fields. When uja has a wire singularity similar to the Dirac monopole in electrody- 
namics, i.e., when the horizon has NUT charge, then Jq would be non-zero.) 

• The mass dipole vanishes, Afi = 0. This tells us that we are in the rest frame of the horizon. 

• For a type I (i.e., spherically symmetric) isolated horizon, Im[4'2] = and Re[^'2] is constant. 
This implies that Mq is the only non-zero multipole moment. 

• If the horizon is symmetric under reflections as in the Kerr solution (i.e., ^2 when 
C 1-^- — C)i then M„ vanishes for odd n while J„ vanishes for even n. 

There is a one-one correspondence between the multipole moments { J„, M„} and the geometry 
of the horizon: Given the horizon area oa and multipoles {J„,M„}, assuming the multipoles 
satisfy a convergence condition for large n, we can reconstruct a non-extremal isolated horizon 
geometry (A,^^^,!?^), uniquely up to diffeomorphisms, such that the area of A is aa and its 
multipole moments are the given {J„,M„}. In vacuum, stationary space-times, the multipole 
moments also suffice to determine the space-time geometry in the vicinity of the horizon. Thus, 
we see that the horizon multipole moments have the expected properties. In the extremal case, 
because of a surprising uniqueness result |148) . the Mn, J„ are universal - the same as those on the 
extremal Kerr IH and the 'true multipoles' which can distinguish one extremal IH from another are 
constructed using different fields in place of ^'2 24 . Finally, note that there is no a-priori reason 
for these source multipoles to agree with the field multipoles at infinity; there could be matter 
fields or radiation outside the horizon which contribute to the field multipoles at infinity. The two 
sets of quantities need not agree even for stationary, vacuum space-times because of contributions 
from the gravitational field in the exterior region. For the Kerr space-time, the source and field 
moments are indeed different for n > 2. However, the difference is small for low n |24) . 

See for further discussion and for the inclusion of electromagnetic fields, and for the 
numerical implementation of these results. 

5.4 Waveform extraction 

The prospect of receiving gravitational waves offers exciting opportunities for astrophysics and for 
testing the dynamical, strong field regime of general relativity (see, e.g., [ZHj)- Black holes are 
among the most promising sources both for terrestrial and space-based observatories. Extracting 
gravitational waveform resulting from, say, gravitational collapse |153| or black hole mergers |138| 
is one of the important goals of numerical relativity. The eventual aim is to provide waveforms 
which, in conjunction with gravitational wave detectors, can be used to study the astrophysics of 
these gravitational wave sources. In this section, we summarize a method of waveform extraction 
based on the isolated horizon formalism. 

The theory of gravitational radiation in exact general relativity is based on structures defined at 
future null infinity 1+ . In particular, associated with every cross-section of I"*" - which represents 
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a retarded 'instant of time' - there is a well defined notion of mass, introduced by Bondi 54) . which 
decreases as gravitational radiation flows across X+ . On the other hand, except for those based 
on conformal methods, most simulations only deal with a finite portion of space-time and thus 
have no direct access to 1+. Instead, one usually uses scalars such as the Weyl tensor component 
^'4 to define the radiation waveform. However, 'i'^ depends on the choice of a null tetrad as 
well as coordinates. While a natural tetrad is available for the perturbation theory on Kerr back 
ground, this is not true in general. In this section we will sketch an approach to solve both these 
problems using the isolated horizon framework: One can construct an approximate analog of X+ 
for a suitable, finite portion of space-time, and introduce a geometrically defined null tetrad and 
coordinates to extract gauge invariant, radiative information from simulations. 




Figure 8: Bondi-like coordinates in a neighborhood of A. 

Our procedure is analogous to the one used at null infinity to construct the Bondi-type coor- 
dinates starting from X+. Let us assume that (A, [£]) is an IH and consider its preferred foliation 
(i.e., the "rest frame" of the horizon) defined in Section [2. 1.31 Using terminology from null infinity, 
we will refer to the leaves of the foliation as the good cuts of A. Pick a null normal from the [£'^]; 
this can be done, e.g., by choosing a specific value for surface gravity. Let Ua be the unique 1-form 
satisfying i^'Ua = — 1 and which is orthogonal to the good cuts. Let {v, 9, (p) be coordinates on A 
such that V is an affine parameter along and the good cuts are given by surfaces of constant v. 
Here and (j> are coordinates on the good cut satisfying Ce9 — and Cecf) = 0. Next, consider past 
null geodesies emanating from the good cuts with —n" as their initial tangent vector. Let the null 
geodesies be affinely parameterized and let the affine parameter be called r and set r = rg on A. 
Lie drag the coordinates {v, 9, (f) along null geodesies. This leads to a set of coordinates (r, u, 6*, 0) 
in a neighborhood of the horizon. The only arbitrariness in this coordinate system is in the choice 
of {9, (j)) on one good cut and the choice of a number tq. Using vacuum Einstein's equations, one 
can obtain a systematic expansion of the metric in inverse powers of — rp) |12| . 

One can also define a null tetrad in the neighborhood in a similar fashion. Let to° be an 
arbitrary complex vector tangent to the good cuts such that (£, n, m, m) is a null tetrad on A. 
Using parallel transport along the null geodesies, we can define a null tetrad in the neighborhood. 
This tetrad is unique up to the spin rotations of m: m — s- e^'^m on the fiducial good cut. This 
construction is shown in Figure |H1 The domain of validity of the coordinate system and tetrads is 
the space-time region in which the null geodesies emanating from good cuts do not cross. 
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Numerically, it might be possible to adapt existing event horizon finders to locate the outgoing 
past null cone of a cross-section of the horizon. This is because event horizon trackers also track 
null surfaces backward in time |82l 183) : the event horizon is the ingoing past null surface starting 
from sufficiently close to future timelike infinity while here we are interested in the outgoing past 
directed null surface starting from an apparent horizon. 

In the Schwarzschild solution, it can be shown analytically that this coordinate system covers 
an entire asymptotic region. In the Kerr space-time, the domain of validity is not explicitly 
known, but in a numerical implementation, this procedure does not encounter geodesic crossing in 
the region of interest to the simulation In general space-times, the extraction of wave forms 
requires the construction to go through only along the past light cones of good cuts which lie in the 
distant future, whence problems with geodesic-crossing are unlikely to prevent one from covering a 
sufficiently large region with these coordinates. This invariantly defined structure provides a new 
approach to extract waveforms. First, the null tetrad presented above can be used to calculate 5*4. 
There is only a phase factor ambiguity (which is a function independent of v and r) inherited from 
the ambiguity in the choice of m° on the fiducial good cut. Second, the past null cone of a good 
cut at a sufficiently late time can be used as an approximate null infinity. This should enable one 
to calculate dynamical quantities such as the analogs of the Bondi mass and the rate of energy loss 
from the black hole, now on the 'approximate' I"*". However, a detailed framework to extract the 
approximate expressions for fluxes of energy and Bondi mass, with sufficient control on the errors, 
is yet to be developed. This is a very interesting analytical problem since its solution would provide 
numerical relativists with an algorithm to extract waveforms and fluxes of energy in gravitational 
waves in an invariant and physically reliable manner. Finally, the invariant coordinates and tetrads 
also enable one to compare late time results of distinct numerical simulations. 
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6 Applications in Mathematical Physics 

For simplicity, in Section 0] we restricted our review of the laws of black hole mechanics to the 
Einstein-Maxwell theory. However, there is a large body of literature on black holes in more 
general theories with dilatonic, Yang-Mills, Higgs, Proca, and Skyrme fields. These fields are not 
expected to be physically significant in the macroscopic, astrophysical world. However, they are of 
considerable interest from a mathematical physics perspective because their inclusion brings about 
qualitative, structural changes in the theory. The most dramatic of these is that the uniqueness 
theorems that play a central role in the Einstein-Maxwell theory are no longer valid. Consequently, 
the structure of these 'colored' or 'hairy' black-holes is much more complicated than those in 
the Einstein-Maxwell theory, and most of the work in this area has been carried out through a 
combination of analytic and numerical methods. By now a very large body of facts about stationary 
black holes with hair has accumulated. A major challenge is to unify this knowledge through a 
few, general principles. 

The isolated horizon framework has provided surprising insights into the properties of hairy 
black holes in equilibrium |19l 1771 1751 1^6 21 76 . While the zeroth and first laws go through in 
a straightforward manner, the notion of the horizon mass now becomes much more subtle and its 
properties have interesting consequences. The framework also suggests a new phenomcnological 
model of colored black holes as bound states of ordinary, uncolored black holes and solitons. This 
model successfully explains the qualitative behavior of these black holes, including their stability 
and instability, and provides unexpected quantitative relations between colored black holes and 
their solitonic analogs. 

In these theories, matter fields are minimally coupled to gravity. If one allows non-minimal 
couplings, the first law itself is modified in a striking fashion: Entropy is no longer given by the 
horizon area but depends also on the matter fields. For globally stationary space-times admitting 
bifurcate Killing horizons, this result was first established by Jacobson, Kang, and Myers |123) . and 
by Iyer and Wald |183II184] for a general class of theories. For scalar fields non-minimally coupled 
to gravity, it has been generalized in the setting of Type II WIHs 22 . While the procedure 
does involve certain technical subtleties, the overall strategy is identical to that summarized in 
Section [4. II Therefore we will not review this issue in detail. 

This section is divided into two parts. In the first, we discuss the mechanics of weakly isolated 
horizons in presence of dilatons and Yang-Mills fields. In the second, we discuss three applications. 
This entire discussion is in the framework of isolated horizons because the effects of these fields on 
black hole dynamics remain largely unexplored. 

6.1 Beyond Einstein— Maxwell theory 

Our primary purpose in this section is to illustrate the differences from the Einstein-Maxwell 
theory. These stem from 'internal charges' and other 'quantum numbers' that are unrelated to 
angular momentum. Therefore, for simplicity, we will restrict ourselves to non-rotating weakly 
isolated horizons. Extension to include angular momentum is rather straightforward. 

6.1.1 Dilatonic couplings 

In dilaton gravity, the Einstein-Maxwell theory is supplemented with a scalar field - called the 
dilaton - and (in the Einstein frame) the Maxwell part of the action is replaced by 



where a > is a free parameter which governs the strength of the coupling of to the Maxwell 
field Fab- If a = 0, one recovers the standard Einstein-Maxwell-Klein-Gordon system, while a = 1 




(69) 
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occurs in a low energy limit of string theory. For our illustrative purposes, it will suffice to consider 
the a — 1 case. 

At spatial infinity, one now has three charges: the ADM mass Madm, the usual electric charge 
Qoo, and another charge Qoo- 

Q. = ^£^, Q^ = lj^e--^^Y. (70) 

Qoo is conserved in space-time (i.e., its value does not change if the 2-spherc of integration is 
deformed) while is not. From the perspective of weakly isolated horizons, it is more useful to 
use OA, Qa, and Qa as the basic charges: 

Qa = ^/;F, QA = ^£e-^^F, (71) 

where S is any cross-section of A. Although the standard electric charge is not conserved in 
space-time, it is conserved along A whence Qa is well-defined. 

It is straightforward to extend the Hamiltonian framework of Section f4.1l to include the dilaton. 
To define energy, one can again seek live time-translation vector fields evolution along which 
is Hamiltonian. The necessary and sufficient condition now becomes the following: There should 
exist a phase space function £'^, constructed from horizon fields, such that 

= ^^"A + ^(t) SQa- (72) 

Thus, the only difference from the Einstein-Maxwell case is that Qa is now replaced by Qa- Again, 
there exists an infinite number of such live vector fields and one can construct them systemati- 
cally starting with any (suitably regular) function /to(flA, Qa) of the horizon area and charge and 
requiring should equal kq. 

The major difference arises in the next step, when one attempts to construct a preferred t^. 
With the dilatonic coupling, the theory has a unique three parameter family of static solutions 
which can be labelled by (oajQajQa) jlOll 1951 I148j . As in the Reissner Nordstrom family, 
these solutions are spherically symmetric. In terms of these parameters, the surface gravity of 
the static Killing field which is unit at infinity, is given by 

The problem in the construction of the preferred is that we need a function kq which depends 
only on oa and Qa and, since depends on all three horizon parameters, one can no longer 
set Kq = on the entire phase space. Thus, there is no live vector field which can generate a 
Hamiltonian evolution and agree with the time-translation Killing field in all static solutions. It 
was the availability of such live vector fields that provided a canonical notion of the horizon mass 
in the Einstein-Maxwell theory in Section Pi. 1.31 

One can weaken the requirements by working on sectors of phase space with fixed values of 
Qa- On each sector, trivially depends only on oa and Qa- So one can set kq = select a 
canonical to, and obtain a mass function Ma = E^a- However, now the first law (|72|) is satisfied 
only if the variation is restricted such that 5Qa = 0. For general variation, one has the modified 
law 1211 

SMa - ^SaA + $ SQa, (74) 

where k = K(to), = {QaQa/Ra) ^'^id Qa — QaQa- Thus, although there is still a first law 
in terms of and Ma, it does not have the canonical form (|7^ because ig does not generate 
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Hamiltonian evolution on the entire phase space. More generally, in theories with multiple scalar 
fields, if one focuses only on static sectors, one obtains similar 'non-standard' forms of the first 
law with new 'work terms' involving scalar fields |1(J0) . From the restricted perspective of static 
sector, this is just a fact. The isolated horizon framework provides a deeper underlying reason: In 
these theories, there is no evolution vector field i° defined for all points of the phase space, which 
coincides with the properly normalized Killing field on all static solutions, and evolution along 
which is Hamiltonian on the full phase space. 

6.1.2 Yang Mills fields 

In the Einstein-Maxwell theory, with and without the dilaton, one can not construct a quantity 
with the dimensions of mass from the fundamental constants in the theory. The situation is different 
for Einstein- Yang-Mills theory because the coupling constant g has dimensions (LA/)^^/^. The 
existence of such a dimcnsionful quantity has interesting consequences. 

ADM Mass 
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Figure 9: The ADM mass as a function of the horizon radius Ra of static spherically symmetric 
solutions to the Einstein-Yang-Mills system (in units provided by the Yang-Mills coupling con- 
stant). Numerical plots for the colorless (n = 0) and families of colored black holes (n — 1,2) are 
shown. (Note that the y-axis begins at M = 0.7 rather than at M — Q.) 



For simplicity, we will restrict ourselves to SU(2) Yang-Mills fields, but results based on the 
isolated horizon framework go through for general compact groups |26) . Let us begin with a sum- 
mary of the known static solutions. First, the Reissner-Nordstrom family constitutes a continuous 
2-parameter set of static solutions of the Einstein- Yang-Mills theory, labelled by (oajQX'^)- 
addition, there is a 1-parameter family of 'embedded Abelian solutions' with (a fixed) magnetic 
charge P^, labelled by (caj-Pa)- Finally, there are families of 'genuinely non-Abelian solutions'. 
For these, the analog of the Israel theorem for Einstein-Maxwell theory fails to hold |127lll28lll29| : 
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the theory admits static solutions which need not be spherically symmetric. In particular, an in- 
finite family of solutions labelled by two integers (rii, 712) is known to exist. All static, spherically 
symmetric solutions are known and they correspond to the infinite sub-family {ni,n2 = 0), la- 
belled by a single integer. However, the two parameter family is obtained using a specific ansatz, 
and other static solutions also exist. Although the available information on the static sector is 
quite rich, in contrast to the Einstein-Maxwell-dilaton system, one is still rather far from having 
complete control. 

However, the existing results are already sufficient to show that, in contrast to the situation in 
the Einstein-Maxwell theory, the ADM mass is not a good measure of the black hole (or horizon) 
mass even in the static case. Let us consider the simplest case, the spherically symmetric static 
solutions labelled by a single integer n (see Figure EJ. Let us decrease the horizon area along any 
branch n ^ 0. In the zero area limit, the solution is known to converge point-wise to a regular, 
static, spherical solution, representing an Einstein- Yang-Mills soliton |38l 11711 16U) . This solution 
has, of course, a non-zero ADM mass M^^M' which equals the limiting value of M^^y^. However, 
in this limit, there is no black hole at alll Hence, this limiting value of the ADM mass can not 
be meaningfully identified with any horizon mass. By continuity, then, M^^j^ can not be taken 
as an accurate measure of the horizon mass for any black hole along an n 7^ branch. Using the 
isolated horizon framework, it is possible to introduce a meaningful definition of the horizon mass 
on any given static branch. 

To establish laws of black hole mechanics, one begins with appropriate boundary conditions. 
In the Maxwell case, the gauge freedom in the vector potential is restricted on the horizon by 
requiring CiAa = on A. The analogous condition ensuring that the Yang-Mills potential A is 
in an 'adapted gauge' on A is more subtle !26;. However, it does exist and again ensures that (i) 
the action principle is well defined, and (ii) the Yang-Mills electric potential $(^) := —\£- A| is 
constant on the horizon, where the absolute sign stands for the norm in the internal space. The 
rest of the boundary conditions are the same as in Section 12.1.11 The proof of the zeroth law 
and the construction of the phase space is now straightforward. There is a well-defined notion of 
conserved horizon charges 

Qr := £ \*F\ d^V, Pl^ :^ -i- £|F| d'V, (75) 

where |F| := [(e''''F^t)(e"''F^j,) with K^j being the Cartan-KiUing metric on SU(2), e"'' the 

alternating tensor on the cross-section S of A, and where |*F| is defined by replacing F with 
*F. Finally, one can again introduce live vector fields t"", time evolution along which generates a 
Hamiltonian flow on the phase space, and establish a first law for each of these t"": 

SEi = ^«(t)^«A + ^it)5QV^. (76) 

Note that, even though the magnetic charge is in general non-zero, it does not enter the 

statement of the first law. In the Abelian case, a non-zero magnetic charge requires non-trivial U(l) 
bundles, and Chcrn numbers characterizing these bundles are discrete. Hence the magnetic charge 
is quantized and, if the phase space is constructed from connections, SP^^ vanishes identically for 
any variation 6. In the non-Abelian case, one can work with a trivial bundle and have non-zero 
pYM^ Therefore, SP^^^ does not automatically vanish and absence of this term in the first law is 
somewhat surprising. 

A more significant difference from the Abelian case is that, because the uniqueness theorem 
fails, one can not use the static solutions to introduce a canonical function kq on the entire phase 
space, whence as in the dilatonic case, there is no longer a canonical horizon mass Ma function 
on the entire phase space. In the next Section 16.21 we will see that it is nonetheless possible to 
introduce an extremely useful notion of the horizon mass for each static sequence. 
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6.2 Structure of colored, static black holes 



We will briefly summarize research in three areas in which the isolated horizon framework has been 
used to illuminate the structure of static, colored black holes and associated solitons''. 

6.2.1 Horizon mass 

Let us begin with Einstein- Yang-Mills theory considered in the last Section I^TI As we saw, the 
ADM mass fails to be a good measure of the horizon mass for colored black holes. The failure 
of black hole uniqueness theorems also prevents the isolated horizon framework from providing a 
canonical notion of horizon mass on the full phase space. However, one can repeat the strategy 
used for dilatonic black holes to define horizon mass unambiguously for the static solutions 77; 

Consider a connected component of the known static solutions, labelled by no = (nj*, nj). Using 
for kq the surface gravity of the properly normalized static Killing vector, in this sector one can 
construct a live vector field and obtain a first law. The energy E"^ is well-defined on the full 
phase space and can be naturally interpreted as the horizon mass Af^"'' for colored black holes 
with 'quantum number' no. The explicit expression is given by 

M'i"^ = —j^pnAR)dR, (77) 

where, following a convention in the literature on hairy black holes |77) . we have used 

(3n„{R) = 2[Kn„{R)]R (78) 

rather than the surface gravity Kf{g{R) of static solutions. Now, as in the Einstein-Maxwell case, 
the Hamiltonian generating evolution along contains only surface terms: Ht^ = i?ADM ~ -^A 
and is constant on each connected, static sector if coincides with the static Killing field on that 
sector. By construction, our ijj has this property for the ng-sector under consideration. Now, 
in the Einstein-Maxwell case, since there is no constant with the dimension of energy, it follows 
that the restriction of ^adm ^'^ static sector must vanish. The situation is quite different in 
Einstein- Yang-Mills theory where the Yang-Mills coupling constant g provides a scale. In c = 1 
units, {gVG)~^ mass. Therefore, we can only conclude that 

^ADM = M^a"^ + (3\/G)-^C("«) (79) 

for some no-dependent constant C'""-'. As the horizon radius shrinks to zero, the static solu- 
tion |127lll8l) under consideration tends to the solitonic solution with the same 'quantum numbers' 
no. Hence, by taking this hmit, we conclude {gVG)~^C^"''^ = M^^^°\ Therefore, on any n-static 
sector, we have the following interesting relation between the black hole and solitonic solutions: 

Ml^il^' = Pn{R) dR + M-^i^^ ^ M^«^('^) + aC'^.) (80) 

Thus, although the main motivation behind the isolated horizon framework was to go beyond 
globally time-independent situations, it has led to an interesting new relation between the ADM 
masses of black holes and their solitonic analogs already in the static sector. 

The relation (|80|) was first proposed for spherical horizons in [77], verified in and extended 
to distorted horizons in |26| . It provided impetus for new work by mathematical physicists working 
on colored black holes. The relation has been confirmed in three more general and non-trivial cases: 

^We are grateful to Alejandro Corichi for correspondence on the recent results in this area. 
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• Non-spherical, non-rotating black holes parameterized by two quantum numbers |1H2| . 

• Non-spherical solutions to the more general Einstein- Yang-Mills-Higgs theory |ll)9) . where 
distortions are caused by 'magnetic dipole hair' |13(J| . 

• Static solutions in the Born-Infeld theory |62|. 

6.2.2 Phenomenological model of colored black holes 

Isolated horizon considerations suggested the following simple heuristic model of colored black 
holes |21| : A colored black hole with quantum numbers n should be thought of as 'bound states' of a 
ordinary (colorless) black hole and a soliton with color quantum numbers n, where ft can be more 
general than considered so far. Thus the idea is that an uncolored black hole is 'bare' and becomes 
'colored' when 'dressed' by the soliton. 

The mass formula (jSOII now suggests that the total ADM mass m'^-j^ has three components: 
the mass A/^""* of the bare horizon, the mass M^^'j^ of the colored soliton, and a binding energy 
given by 

^binding = Mi"'(i?A) - M^^\Ra). (81) 

If this picture is correct, being gravitational binding energy, -Ebinding would have to be negative. 
This expectation is borne out in explicit examples. The model has several predictions on the qual- 
itative behavior of the horizon mass (|77|) . the surface gravity, and the relation between properties 
of black holes and solitons. We will illustrate these with just three examples (for a complete list 
with technical caveats, see |21)'l: 

• For any fixed value of i?A and of all quantum numbers except n, the horizon mass and surface 
gravity decrease monotonically with n. 

• For fixed values of all quantum numbers n, the horizon mass M^\Ra) is non-negative, 
vanishing if and only if i?A vanishes, and increases monotonically with i?A- /?(«) is positive 
and bounded above by 1. 

• For fixed n, the binding energy decreases as the horizon area increases. 

The predictions for fixed n have recently been verified beyond spherical symmetry: for the 
distorted, axially symmetric Einstein- Yang-Mills solutions in |132| and for the distorted 'dipole' 
solutions in Einstein- Yang-Mills-Higgs solutions in jlU9) . Taken together, the predictions of this 
model can account for all the qualitative features of the plots of the horizon mass and surface 
gravity as functions of the horizon radius and quantum numbers. More importantly, they have 
interesting implications on the stability properties of colored black holes. 

One begins with an observation about solitons and deduces properties of black holes. Einstein- 
Yang-Mills solitons are known to be unstable (173) : under small perturbations, the energy stored 
in the 'bound state' represented by the soliton is radiated away to future null infinity X+. The 
phenomenological model suggests that colored black holes should also be unstable and they should 
decay into ordinary black holes, the excess energy being radiated away to infinity. In general, 
however, even if one component of a bound system is unstable, the total system may still be stable 
if the binding energy is sufficiently large. An example is provided by the deuteron. However, an 
examination of energetics reveals that this is not the case for colored black holes, so instability 
should prevail. Furthermore, one can make a few predictions about the nature of instability. We 
summarize these for the simplest case of spherically symmetric, static black holes for which there 
is a single quantum number n: 
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1 . All the colored black holes on the nth branch have the same number (namely, 2n) of unstable 
modes as the nth soliton. (The detailed features of these unstable modes can differ especially 
because they are subject to different boundary conditions in the two cases.) 

2. For a given n, colored black holes with larger horizon area are less unstable. For a given 
horizon area, colored black holes with higher value of n are more unstable. 

3. The 'available energy' for the process is given by 

E^^^.^-M^-\El^lSL^- (82) 

Part of it is absorbed by the black hole so that its horizon area increases and the rest is 
radiated away to infinity. Note that -E^vaiiabie ^® computed knowing just the initial 
configuration. 

4. In the process the horizon area necessarily increases. Therefore, the energy radiated to 
infinity is strictly less than ^^^"iii^bic- 

Expectation ^ of the model is known to be correct jl72| . Prediction El has been shown to be 
correct in the n = 1, colored black holes in the sense that the frequency of all unstable modes is 
a decreasing function of the area, whence the characteristic decay time grows with area (181. .51] . 
To our knowledge a detailed analysis of instability, needed to test Predictions 13 and 0] are yet to 
be made. 

Finally, the notion of horizon mass and the associated stability analysis has also provided an 
'explanation' of the following fact which, at first sight, seems puzzling. Consider the 'embedded 
Abelian black holes' which are solutions to Einstein- Yang-Mills equations with a specific magnetic 
charge P^. They are isometric to a family of magnetically charged Reissner-Nordstrom solutions 
and the isometry maps the Maxwell field strength to the Yang-Mills field strength. The only 
difference is in the form of the connection; while the Yang-Mills potential is supported on a trivial 
SU(2) bundle, the Maxwell potential requires a non-trivial U(l) bundle. Therefore, it comes as 
an initial surprise that the solution is stable in the Einstein-Maxwell theory but unstable in the 
Einstein- Yang-Mills theory |52II61| . It turns out that this difference is naturally explained by the 
WIH framework. Since the solutions are isometric, their ADM mass is the same. However, since 
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the horizon mass arises from Hamiltonian considerations, it is theory dependent: It is lower in the 
Einstein- Yang-Mills theory than in the Einstein-Maxwell theory! Thus, from the Einstein- Yang- 
Mills perspective, part of the ADM mass is carried by the soliton and there is positive -E^aUabic 
which can be radiated away to infinity. In the Einstein-Maxwell theory, i?^^iiabio zero. The 
stability analysis sketched above therefore implies that the solution should be unstable in the 
Einstein- Yang-Mills theory but stable in the Einstein-Maxwell theory. This is another striking 
example of the usefulness of the notion of the horizon mass. 

6.2.3 More general theories 

We will now briefly summarize the most interesting result obtained from this framework in more 
general theories. When one allows Higgs or Proca fields in addition to Yang-Mills, or considers 
Einstein-Skyrme theories, one acquires additional dimensionful constants which trigger new phe- 
nomena | 5U[ 11781 1181) . One of the most interesting is the 'crossing phenomena' of Figure [TTI where 
curves in the 'phase diagram' (i.e., a plot of the ADM mass versus horizon radius) corresponding 
to the two distinct static families cross. This typically occurs in theories in which there is a length 
scale even in absence of gravity, i.e., even when Newton's constant is set equal to zero [1541 121| . 




Figure 11: The ADM mass as a function of the horizon radius Ra in theories with a built-in 
non- gravitational length scale. The schematic plot shows crossing of families labelled by n = 1 and 
n = 2 at Ra^ R^a°'. 

In this case, the notion of the horizon mass acquires further subtleties. If, as in the Einstein- 
Yang-Mills theory considered earlier, families of static solutions carrying distinct quantum numbers 
do not cross, there is a well-defined notion of horizon mass for each static solution, although, as the 
example of 'embedded Abelian solutions' shows, in general its value is theory dependent. When 
families cross, one can repeat the previous strategy and use Equation (|77|) to define a mass 
along each branch. However, at the intersection point i?^*''"' of the nth and (n -I- l)th branches, 
the mass is discontinuous. This discontinuity has an interesting implication. Consider the closed 
curve C in the phase diagram, starting at the intersection point and moving along the (n + l)th 
branch in the direction of decreasing area until the area becomes zero, then moving along oa = 
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to the nth branch and moving up to the intersection point along the nth branch (see Fierure [TT|l . 
Discontinuity in the horizon mass implies that the integral of /3(i?A) along this closed curve is 
non-zero. Furthermore, the relation between the horizon and the soliton mass along each branch 
implies that the value of this integral has a direct physical interpretation: 

MiT'^-Mi:^ = ^fjir)dr. (83) 

This is a striking prediction because it relates differences in masses of solitons to the knowledge 
of horizon properties of the corresponding black holesl Because of certain continuity properties 
which hold as one approaches the static Einstein- Yang-Mills sector in the space to static solutions 
to Einstein- Yang-Mills-Higgs equations |7B], one can also obtain a new formula for the ADM 
masses of Einstein- Yang-Mills solitons. If (1 — /3(„-)(r)) for the black hole solutions of this theory 
is integrable over the entire positive half line, one has |21) : 

M^:i = ^j^{l-PUr))dr. (84) 

Both these predictions of the phenomenological model [2J have been verified numerically in the 
spherically symmetric case [76) . but the axi-symmetric case is still open. 
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7 Applications in Quantum Gravity 



As discussed in the introduction, laws of black hole mechanics, discovered in the early seventies, 
provided a concrete challenge to candidate quantum theories of gravity: Account for the ther- 
modynamic, black hole entropy through a detailed, statistical mechanical counting of appropriate 
micro-states. Indeed, this is essentially the only concrete quantitative hint we have had about 
the nature of quantum space-time geometry. The isolated horizon framework has been used to 
address this issue systematically and has led to the only available detailed calculations within a 
full-fledged approach to quantum gravity that encompass realistic black holes (which carry no or 
negligible gauge charges and may be distorted) . As we will discuss in Section |S1 what we know 
about dynamical horizons does suggest that there should be interesting generalizations of these 
results also to non-equilibrium situations. But so far there has been no work along these lines. 

7.1 Preliminaries 

The first and the second laws of black hole mechanics discussed in Section 01 

5E = -^(5a + work, d6a > 0, (85) 

have a close similarity with the first and second laws of thermodynamics 

SE = T5S + work, 6S > 0. (86) 

This suggest that one should assign to black holes an entropy proportional to their area. Indeed, 
already in the early seventies, using imaginative thought experiments, Bekenstein |44l I45| argued 
that such an assignment is forced upon us if, in presence of black holes, we are not going to simply 
abandon the second law of thermodynamics. Specifically, consider an experiment in which a box 
containing radiation is slowly lowered into a black hole. At the end of the process, the entropy of the 
exterior region decreases in apparent violation of the second law. Bekenstein argued that the area 
of the horizon increases in the process so that the total entropy would not decrease if the black hole 
is assigned entropy equal to an appropriate multiple of the area. Hawking's celebrated calculation 
of black hole evaporation provided the precise numerical coefficient. For, it suggested that 
we should assign the black hole a temperature T = hn/i'K. By comparing the forms of the first 
laws of black hole mechanics and thermodynamics, one is then led to set 

= with 4i = Gh. (87) 

Note that the factor of h is absolutely essential: In the limit h ^ the black hole temperature 
goes to zero and its entropy tends to infinity just as one would expect classically. 

The relation 1)87(1 is striking and deep because it brings together the three pillars of fundamental 
physics - general relativity, quantum theory, and statistical mechanics. However, the argument 
itself is a rather hodge-podge mixture of classical and semi-classical ideas, reminiscent of the Bohr 
theory of atom. A natural question then is: What is the analog of the more fundamental, Pauli- 
Schrodinger theory of the hydrogen atom? More precisely, what is the statistical mechanical origin 
of black hole entropy? To answer this question, one has to isolate the microscopic degrees of 
freedom responsible for this entropy. For a classical ideal gas, these are given by the positions 
and momenta of all molecules; for a magnet, by the states of each individual spin at lattice sites. 
How do we represent the analogous microscopic degrees of freedom for a black hole? They can not 
be described in terms of quantum states of physical gravitons because we are dealing with black 
holes in equilibrium. In the approach based on weakly isolated horizons, they are captured in 
the quantum states of the horizon geometry. Just as one would expect from Bekenstein's thought 
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experiments, these degrees of freedom can interact with the exterior curved geometry, and the 
resulting notion of black hole entropy is tied to observers in the exterior regions. 

A heuristic framework for the calculation of entropy was suggested by John Wheeler in the 
nineties, which he christened ^It from bit\ Divide the black hole horizon into elementary cells, 
each with one Planck unit of area -^pj and assign to each cell two microstates. Then the total 
number of states Af is given hy Af = 2", where n = (a/£p[) is the number of elementary cells, 
whence entropy is given by S* = \nj\f a. Thus, apart from a numerical coefficient, the entropy 
{'IV) is accounted for by assigning two states to each elementary cell. This qualitative 

picture is simple and attractive. But it raises at least three central questions: 

• What is the origin of the 'elementary cells' and why is each endowed with an area ip-^? 

• What is the origin of the microstates carried by each elementary cell and why are there 
precisely two microstates? 

• What does all this have to do with a black hole? Why doesn't it apply to any 2-surface, 
including a 2-sphere in Minkowski space-time? 

An understanding of geometry of quantum WIHs provides a detailed framework which, in 
particular, answers these questions. The precise picture, as usual, is much more involved than that 
envisaged by Wheeler. Eigenvalues of area operator turn out to be discrete in quantum geometry 
and one can justify dividing the horizon 2-sphere into elementary cells. However, there are many 
permissible area eigenvalues and cells need not all carry the same area. To identify horizon surface 
states that are responsible for entropy, one has to crucially use the WIH boundary conditions. 
However, the number of surface states assigned to each cell is not restricted to two. Nonetheless, 
Wheeler's picture turns out to be qualitatively correct. 

7.2 Quantum horizon geometry 

For simplicity of presentation, in this section we will restrict ourselves to type I WIHs, i.e., the ones 
for which the only non-zero multipole moment is the mass monopole. The extension to include 
type H horizons with rotations and distortion will be briefly summarized in the next Section 17.31 
Details can be found in [TUl [TTl QUI 151 El ] . 

The point of departure is the classical Hamiltonian formulation for space-times A4 with a type 
I WIH A as an internal boundary, with flxed area oq and charges Qq, where a runs over the 
number of distinct charges (Maxwell, Yang-Mills, dilaton, . . . ) allowed in the theory. As we 
noted in Section 14.11 the phase space T can be constructed in a number of ways, which lead to 
equivalent Hamiltonian frameworks and first laws. However, so far, the only known way to carry 
out a background independent, non-perturbative quantization is through connection variables 

As in Figure let us begin with a partial Cauchy surface M whose internal boundary in A4 
is a 2-sphere cross-section S* of A and whose asymptotic boundary is a 2-sphere 6*00 at spatial 
infinity. The configuration variable is an SU(2) connection on M, where i takes values in the 
3-dimensional Lie-algebra su(2) of SU(2). Just as the standard derivative operator acts on tensor 
fields and enables one to parallel transport vectors, the derivative operator constructed from 
acts on fields with internal indices and enables one to parallel transport spinors. The conjugate 
momentum is represented by a vector field Pf with density weight 1 which also takes values in 
su(2); it is the analog of the Yang-Mills electric field. (In absence of a background metric, momenta 
always carry a density weight 1.) can be regarded as a (density weighted) triad or a 'square- root' 
of the intrinsic metric qab on S: {8TTGj)^P°'Pj5ij = qq"'^, where 6ij is the Cartan Killing metric 
on su(2), q is the determinant of ijab and 7 is a positive real number, called the Barbero-Immirzi 
parameter. This parameter arises because there is a freedom in adding to Palatini action a multiple 
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of the term which is 'dual' to the standard one, which does not affect the equations of motion but 
changes the definition of momenta. This multiple is 7. The presence of 7 represents an ambiguity 
in quantization of geometry, analogous to the 0-ambiguity in QCD. Just as the classical Yang- 
Mills theory is insensitive to the value of 9 but the quantum Yang-Mills theory has inequivalent 
0-sectors, classical relativity is insensitive to the value of 7 but the quantum geometries based on 
different values of 7 are (unitarily) inequivalent (for details, see, e.g., |82|). 

Thus, the gravitational part of the phase space F consists of pairs (A^,Pf) of fields on M 
satisfying the boundary conditions discussed above. Had there been no internal boundary, the 
gravitational part of the symplectic structure would have had just the expected volume term: 



J M 



I M 

where Sahi '■= VabcPi is the 2-form dual to the momentum and Si, and 62 denote any two 
tangent vectors to the phase space. However, the presence of the internal boundary changes the 
phase space-structure. The type I WIH conditions imply that the non-trivial information in the 
pull-back Al^ of to S is contained in a U(l) connection Wa ■— Al^r.i on 5*, where is the unit, 
internal, radial vector field on S. Its curvature 2-form Fab is related to the pull-back of the 
2-form S^^ to S via 

F = dW^ ~8nGj— SV,. (89) 
ao ^ 

The restriction is called the horizon boundary condition. In the Schwarzschild space-time, all the 
2-spheres on which it is satisfied lie on the horizon. Finally, the presence of the internal boundary 
modifies the symplectic structure: It now acquires an additional boundary term 

n{Si,S2) ^nv{Si,S2) + ns{du62), (90) 

with 

ns{Si,62)^^£^j^JiWA52W, (91) 

The new surface term is precisely the symplectic structure of a well-known topological field theory 
- the U(l)-Chern-Simons theory. The symplectic structures of the Maxwell, Yang-Mills, scalar, 
and dilatonic fields do not acquire surface terms. Conceptually, this is an important point: This, 
in essence, is the reason why (for minimally coupled matter) the black hole entropy depends just 
on the area and not, in addition, on the matter charges. 

In absence of internal boundaries, the quantum theory has been well- understood since the 
mid-nineties (for recent reviews, see, [1621 11751 132) ). The fundamental quantum excitations are 
represented by Wilson lines (i.e., holonomies) defined by the connection and are thus 1-dimensional, 
whence the resulting quantum geometry is polymer-like. These excitations can be regarded as flux 
lines of area for the following reason. Given any 2-surface S on M , there is a self-adjoint operator 
A§ all of whose eigenvalues are known to be discrete. The simplest eigenvectors are represented by 
a single flux line, carrying a half-integer j as a label, which intersects the surface S exactly once, 
and the corresponding eigenvalue ag of is given by 



as = 87r74iv/j(7 + T). (92) 

Thus, while the general form of the eigenvalues is the same in all 7-sectors of the quantum theory, 
their numerical values do depend on 7. Since the eigenvalues are distinct in different 7-sectors, it 
immediately follows that these sectors provide unitarily inequivalent representations of the algebra 
of geometric operators; there is 'super-selection'. Put differently, there is a quantization ambiguity, 
and which 7-sector is actually realized in Nature is an experimental question. One appropriate 
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experiment, for example a measurement of the smallest non-zero area eigenvalue, would fix the 
value of 7 and hence the quantum theory. Every further experiment - e.g., the measurement of 
higher eigenvalues or eigenvalues of other operators such as those corresponding to the volume of a 
region - would provide tests of the theory. While such direct measurements are not feasible today 
we will see that, somewhat surprisingly, the Hawking-Bekenstein formula H87(l for the entropy of 
large black holes provides a thought experiment to fix the value of 7. 

Recall next that, because of the horizon internal boundary, the symplectic structure now has 
an additional surface term. In the classical theory, since all fields are smooth, values of fields on 
the horizon are completely determined by their values in the bulk. However, a key point about 
field theories is that their quantum states depend on fields which are arbitrarily discontinuous. 
Therefore, in quantum theory, a decoupling occurs between fields in the surface and those in the 
bulk, and independent surface degrees of freedom emerge. These describe the geometry of the 
quantum horizon and are responsible for entropy. 

In quantum theory, then, it is natural to begin with a total Hilbert space Ti. = Tiy <E) Ti-s where 
Tiv is the well-understood bulk or volume Hilbert space with 'polymer-like excitations', and Tis is 
the surface Hilbert space of the U(l)-Chern-Simons theory. As depicted in Figure IT^ the polymer 
excitations puncture the horizon. An excitation carrying a quantum number j 'deposits' on S an 
area equal to Sirip^ \/ j{j + !)■ These contributions add up to endow S a total area oq. The surface 
Chern-Simons theory is therefore defined on the punctured 2-sphere S. To incorporate the fact 
that the internal boundary S is not arbitrary but comes from a WIH, we still need to incorporate 
the residual boundary condition (|89(l . This key condition is taken over as an operator equation. 
Thus, in the quantum theory, neither the triad nor the curvature of W are frozen at the horizon; 
neither is a classical field. Each is allowed to undergo quantum fluctuations, but the quantum 
horizon boundary condition requires that they have to fluctuate in tandem. 



Figure 12: Quantum horizon. Polymer excitations in the bulk puncture the horizon, endowing 
it with quantized area. Intrinsically, the horizon is flat except at punctures where it acquires a 
quantized deficit angle. These angles add up to endow the horizon with a 2-sphere topology. 

An important subtlety arises because the operators corresponding to the two sides of Equa- 
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tion H89|l act on different Hilbert spaces: While F is defined on Hs, S • r is defined on Hv- 
Therefore, the quantum horizon boundary condition introduces a precise intertwining between the 
bulk and the surface states: Only those states 5* = (E) in Ti. which satisfy 

(1(8)/')^'= ^-^ET«)f^ 5- (93) 

can describe quantum geometries with WIH as inner boundaries. This is a stringent restriction: 
Since the operator on the left side acts only on surface states and the one on the right side acts only 
on bulk states, the equation can have solutions only if the two operators have the same eigenvalues 
(in which case we can take to be the tensor product of the two eigenstates). Thus, for solutions 
to Equation H93|l to exist, there has to be a very delicate matching between eigenvalues of the triad 
operators S • r calculated from bulk quantum geometry, and eigenvalues of F, calculated from 
Chern-Simons theory. The precise numerical coefficients in the surface calculation depend on the 
numerical factor in front of the surface term in the symplectic structure, which is itself determined 
by our type I WIH boundary conditions. Thus, the existence of a coherent quantum theory of WIHs 
requires that the three cornerstones - classical isolated horizon framework, quantum mechanics of 
bulk geometry, and quantum Chern-Simons theory - be united harmoniously. Not only should the 
three frameworks be mutually compatible at a conceptual level, but certain numerical coefficients, 
calculated independently within each framework, have to match delicately. Remarkably, these 
delicate constraints are met, whence the quantum boundary conditions do admit a sufficient number 
of solutions. 

We will conclude by summarizing the nature of geometry of the quantum horizon that results. 
Given any state satisfying Equation H93|l . the curvature of vanishes everywhere except at 
the points at which the polymer excitations in the bulk puncture S. The holonomy around each 
puncture is non-trivial. Consequently, the intrinsic geometry of the quantum horizon is flat except 
at the punctures. At each puncture, there is a deficit angle, whose value is determined by the 
holonomy of W around that puncture. Each deficit angle is quantized and these angles add up to 
Att as in a discretized model of a 2-sphere geometry. Thus, the quantum geometry of a WIH is 
quite different from its smooth classical geometry. 

7.3 Entropy 

Let us now summarize the ideas behind counting of surface microstates that leads to the expression 
of entropy. To incorporate dynamics in this canonical approach, we have to first construct physical 
states by imposing quantum Einstein equations (i.e., quantum constraints). While the procedure 
is technically quite involved, the result is simple to state: What matters is only the number of 
punctures and not their locations. To calculate entropy, then, one constructs a micro-canonical 
ensemble as follows. Fix the number n of punctures and allow only those (non-zero) spin-labels ji 
and charge labels qf on the polymer excitations which endow the horizon with a total area in an 
interval (oq — e, oq -t- e) and charges in an interval {Qq — e", Qq + e"). (Here / = 1,2, ... ,n and e 
and e" are suitably small. Their precise values will not affect the leading contribution to entropy.) 
We denote by the sub-space oiH = Hv^'Hs in which the volume states ^'y are chosen with 
the above restrictions on jj and qf , and the total state 4* satisfies the quantum horizon boundary 
condition as well as quantum Einstein equations. Then the desired micro-canonical ensemble 
consists of states in = ©„7iJ^^. Note that, because there is no contribution to the symplectic 
structure from matter terms, surface states in Ti^^ refer only to the gravitational sector. 

The next step is to calculate the entropy of this quantum, micro-canonical ensemble. Note 
first that what matters are only the surface states. For, the 'bulk-part' describes, e.g., states of 
gravitational radiation and matter fields far away from A and are irrelevant for the entropy Sa of 
the WIH. Heuristically, the idea then is to 'trace over' the bulk states, construct a density matrix 



63 



jOBH describing a maximum-entropy mixture of surface states and calculate tr pen In pbh ■ As is 
usual in entropy calculations, this translates to the evaluation of the dimension AT of a well-defined 
sub-space 7i™ of the surface Hilbert space, namely the linear span of those surface states which 
occur in TC^^. Entropy Sa is given by In A/". 

A detailed calculation |84II149| leads to the following expression of entropy: 

5. = - i In ^ + O (in i|) , with ,0 0.2375, (94) 

where O(hi{ao/£p^)) is a term which, when divided by ao/^pj, tends to zero as ao/£p tends to infinity. 
Thus, for large black holes, the leading term in the expression of entropy is proportional to area. 
This is a non-trivial result. For example if, as in the early treatments |17(JI 1361 11631 Ilb4l 11341 1135| 
one ignores the horizon boundary conditions and the resulting Chern-Simons term in the symplectic 
structure, one would find that Sa is proportional to However, the theory does not have a 

unique prediction because the numerical coefficient depends on the value of the Barbero-Immirzi 
parameter 7. The appearance of 7 can be traced back directly to the fact that, in the 7-sector of 
the theory, the area eigenvalues are proportional to 7. 

One adopts a 'phenomenological' viewpoint to fix this ambiguity. In the infinite dimensional 
space of geometries admitting A as their inner boundary, one can fix one space-time, say the 
Schwarzschild space-time with mass Mq ^ Mpi, (or, the de Sitter space-time with the cosmological 
constant Aq <C l/^pj, or, ... ). For agreement with semi-classical considerations in these cases, the 
leading contribution to entropy should be given by the Hawking-Bekenstein formula (|87() . This 
can happen only in the sector 7 = 70- The quantum theory is now completely determined through 
this single constraint. We can go ahead and calculate the entropy of any other type I WIH in this 
theory. The result is again Sa — Sbh- Furthermore, in this 7-sector, the statistical mechanical 
temperature of any type I WIH is given by Hawking's semi-classical value kH/ (27r) |^I135| . Thus, 
we can do one thought experiment - observe the temperature of one large black hole from far away 
- to eliminate the Barbero-Immirzi ambiguity and fix the theory. This theory then predicts the 
correct entropy and temperature for all WIHs with oq ^ irrespective of other parameters such 
as the values of the electric or dilatonic charges or the cosmological constant. An added bonus 
comes from the fact that the isolated horizon framework naturally incorporates not only black hole 
horizons but also the cosmological ones for which thermodynamical considerations are also known 
to apply innj- The quantum entropy calculation is able to handle both these horizons in a single 
stroke, again for the same value 7 = 70 of the Barbero-Immirzi parameter. In this sense, the 
prediction is robust. 

Finally, these results have been subjected to further robustness tests. The first comes from non- 
minimal couplings. Recall from Section that in presence of a scalar field which is non-minimally 
coupled to gravity, the first law is modified [1231 11831 1184) . The modification suggests that the 
Hawking-Bekenstein formula S'bh = oo/(4£p[) is no longer valid. If the non- minimal coupling is 
dictated by the action 

S[gab, <l^]^ j {j^fWR ~ ig'^^V.^Vb^ - y (0)) d^x, (95) 

where R is the scalar curvature of the metric gab and V is a potential for the scalar field, then the 
entropy should be given by |22 

^A = ^^/(0)rf'K (96) 

An immediate question arises: Can the calculation be extended to this qualitatively new situation? 
At first, this seems very difficult within an approach based on quantum geometry, such as the one 
described above, because the relation is non-geometric. However, the answer was shown to be in 
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the affirmative |2()| for type I horizons. It turns out that the metric qab on M is no longer coded 
in the gravitational momentum alone but depends also on the scalar field. This changes the 
surface term fig in the symplectic structure H9U|I as well as the horizon boundary condition (|89|l 
just in the right way for the analysis to go through. Furthermore, state counting now leads to the 
desired expression H96|l precisely for the same value 7 = 70 of the Barbero-Immirzi parameter |2()| . 

Next, one can consider type II horizons which can be distorted and rotating. In this case, 
all the (gravitational, electro-magnetic, and scalar field) multipoles are required as macroscopic 
parameters to fix the system of interest. Therefore, now the appropriate ensemble is determined 
by fixing all these multipoles to lie in a small range around given values. This ensemble can be 
constructed by first introducing multipole moment operators and then restricting the quantum 
states to lie in the subspace of the Hilbert space spanned by their eigenvectors with eigenvalues 
in the given intervals. Again recent work shows that the state counting yields the Hawking- 
Bekenstein formula (|87|) for minimally coupled matter and its modification 19t)|l for non-minimally 
coupled scalar field, for the same value 7 = 70 of the Barbero-Immirzi parameter pil25|. 

To summarize, the isolated horizon framework serves as a natural point of departure for a 
statistical mechanical calculation of black hole entropy based on quantum geometry. How does 
this detailed analysis compare with the 'It from Bit' scenario )187| with which we began? First, the 
quantum horizon boundary conditions play a key role in the construction of a consistent quantum 
theory of the horizon geometry. Thus, unlike in the 'It from Bit' scenario, the calculation pertains 
only to those 2-spheres S which are cross-sections of a WIH. One can indeed divide the horizon 
into elementary cells as envisaged by Wheeler: Each cell contains a single puncture. However, the 
area of these cells is not fixed but is dictated by the j-label at the puncture. Furthermore, there 
are not just 2 but rather 2j + 1 states associated with each cell. Thus, the complete theory is much 
more subtle than that envisaged in the 'It from Bit' scenario. 
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8 Outlook 



In the last six sections, we summarized the isolated and dynamical horizon frameworks and their 
applications. These provide a quasi- local and more physical paradigm for describing black holes 
both in the equilibrium and dynamical regimes. One of the most pleasing aspects of the paradigm 
is that it provides a unified approach to a variety of problems involving black holes, ranging from 
entropy calculations in quantum gravity, to analytical issues related to numerical simulations, to 
properties of hairy black holes, to laws of black hole mechanics. More importantly, as summarized 
in Section^ these frameworks enable one to significantly extend the known results based on Killing 
and event horizons, and provide brand new laws in the dynamical regime. 

In this section, we will discuss some of the open issues whose resolution would lead to significant 
progress in various areas. 

Isolated horizons 

This is the best understood part of the new paradigm. Nonetheless, several important issues 
still remain. We will illustrate these with a few examples: 

Black hole mechanics 

Throughout, we assumed that the space-time metric is C'' (with > 3) and the topology 
of A is 5*^ X R. These assumptions rule out, by fiat, the presence of a NUT charge. 
To incorporate a non-zero NUT charge in black hole mechanics, one must either allow 
A to be topologically S^, or allow for 'wire singularities' in the rotation 1-form tOa on 
A. The zeroth law goes through in these more general situations. What about the first 
law? Arguments based on Euclidean methods |114l 11461 1147| show that entropy is no 
longer given by the horizon area but there is also contribution due to 'Misner strings'. 
However, to our knowledge, a systematic derivation in the Lorentzian regime is not yet 
available. Such a derivation would provide a better understanding of the physical origin 
of the extra terms. The covariant phase space methods used in the isolated horizon 
framework should be applicable in this case. 
Application to numerical relativity 

We saw in Section [5.41 that, in the IH framework, one can introduce an approximate 
analog of future null infinity and invariant coordinate systems and tetrads in its 
neighborhood. With this structure, it is feasible to extract waveforms and energy fluxes 
in a reliable manner within the standard 3-1-1 Cauchy evolution of numerical relativity, 
without having to do a Cauchy characteristic matching or use conformal field equations. 
The challenge here is to develop, on the approximate X"*", the analog of the basics of the 
Bondi framework [MJ 1185| at null infinity. 

Colored black holes 

As discussed in Sectional black hole uniqueness theorems of the Einstein-Maxwell the- 
ory fail once non-Abelian fields are included. For example, there are black hole solutions 
to the Einstein- Yang-Mills equations with non-trivial Yang-Mills fields whose only non- 
zero charge at infinity is the ADM mass. Thus, from the perspective of infinity, they 
are indistinguishable from the Schwarzschild solution. However, their horizon properties 
are quite different from those of the Schwarzschild horizon. This example suggests that 
perhaps the uniqueness theorems fail because of the insistence on evaluating charges at 
infinity. Corichi, Nucamendi, and Sudarsky [TSj have conjectured that the uniqueness 
theorems could be restored if they are formulated in terms of all the relevant horizon 
charges. This is a fascinating idea and it has been pursued numerically. However, care is 
needed because the list of all relevant charges may not be obvious a priori. For example, 
in the case of static but not necessarily spherical Yang-Mills black holes, the conjec- 
ture seemed to fail |132) until one realized that, in addition to the standard Yang-Mills 
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charges at the horizon, one must also include a topological, 'winding charge' in the 
list Once this charge is included, uniqueness is restored not only in the static sector 
of the Einstein- Yang-Mills theory, but also when Higgs fields '109" and dilatons 
are included. The existence of these semi-numerical proofs suggests that it should be 
possible to establish uniqueness completely analytically. A second set of problems in- 
volves the surprising relations, e.g., between the soliton masses and horizon properties 
of colored black holes, obtained using isolated horizons. Extensions of these results to 
situations with non-zero angular momentum should be possible and may well provide 
yet new insights. 

Quantum black holes 

In the approach based on isolated horizons, the microscopic degrees of freedom re- 
sponsible for the statistical mechanical entropy of black holes are directly related to the 
quantum geometry of horizons. Therefore, their relation to the curved space-time geom- 
etry is clearer than in, say, the string theory calculations based on D-branes. Therefore, 
one can now attempt to calculate the Hawking radiation directly in the physical space- 
time as a process in which quanta of area are converted to quanta of matter. However, 
such a calculation is yet to be undertaken. A direct approach requires quantum field 
theory (of matter fields) on a quantum geometry state which is approximated by the 
classical black hole space-time. Elements of this theory are now in place. The concrete 
open problem is the calculation of the absorption cross-section for quantum matter fields 
propagating on this 'background state'. If this can be shown to equal the classical ab- 
sorption cross-section to the leading order, it would follow (22] that the spectrum of the 
outgoing particles would be thermal at the Hawking temperature. Another, perhaps 
more fruitful, avenue is to introduce an effective model whose Hamiltonian captures the 
process by which quanta of horizon area are converted to quanta of matter. Both these 
approaches are geared to large black holes which can be regarded as being in equilibrium 
for the process under consideration, i.e., when oa 2> ^pj. However, this approximation 
would fail in the Planck regime whence the approaches can not address issues related 
to 'information loss'. Using ideas first developed in the context of quantum cosmology, 
effects of the quantum nature of geometry on the black hole singularity have recently 
been analyzed |17| . As in the earlier analysis of the big-bang singularity, it is found that 
the black hole singularity is resolved, but the classical space-time dissolves in the Planck 
regime. Therefore, the familiar Penrose diagrams are no longer faithful representations 
of the physical situation. Suppose that, evolving through what was singularity in the 
classical theory, the quantum state becomes semi-classical again on the 'other side'. 
Then, the indications are that information would not be lost: It would be recovered 
on full , although observers restricted to lie in the part of space-time which is com- 
pletely determined by the data on T~ would see an approximate Hawking radiation. If 
on the other hand the evolved state on the 'other side' never becomes semi-classical, 
information would not be recovered on the available T'^ . An outstanding open issue is 
which of this possibility is actually realized. 

Dynamical horizons 

The DH framework is less developed and the number of open issues is correspondingly higher. 
At least for the classical applications of the framework, these problems are more important. 

Free data and multipoles 

Since H is space-like, to find the fields {qab,Kab) which constitute the DH geometry, 
one has to solve just the initial value equations, subject to the condition that H admits 
a foliation by marginally trapped surfaces. A general solution of this problem would 
provide the 'free data' on DHs. In the case when the marginally trapped surfaces are 
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round spheres, this problem has been analyzed by Bartnik and Isenberg pTT] , As noted 
in Section in this case there are no DHs in absence of matter sources. In presence 
of matter, one can freely specify the trace K of the extrinsic curvature and the (radial 
component of the) momentum density Tabf^f^ , and determine the geometry by solving 
a non-linear ordinary differential equation whose solutions always exist locally. It would 
be interesting to analyze the necessary and sufficient conditions on the free data which 
guarantee that global solutions exist and the DH approaches the Schwarzschild horizon 
asymptotically. From the point of view of numerical relativity, a more pressing chal- 
lenge is to solve the constraint equations in the vacuum case, assuming only that the 
marginally trapped surfaces are axi-symmetric. Using the free data, as in the case of 
isolated horizons one could define multipoles. Since ^"2 is again defined unambigu- 
ously, a natural starting point is to use it as the key geometrical object as in the IH case. 
However, just as the Bondi mass aspect acquires shear terms in presence of gravitational 
radiation [SI], it is likely that, in the transition from isolated to dynamical horizons, ^'2 
would have to be supplemented with terms involving, e.g., aat (and perhaps also C°)- 
For instance, by adding a suitable combination of such terms, one may be able to relate 
the rate of change of the mass quadrupole moment with the flux of energy across H. 
Geometric analysis 

The dynamical horizon framework provides new inputs for the proof of Penrose inequal- 
ities which, when applied to time symmetric data (i.e., when the extrinsic curvature 
vanishes), say that the total (ADM) mass of space-time must be greater than half the 
radius of the apparent horizon on any Cauchy slice. This conjecture was recently proved 
by Bray |59j . and by Huiskcn and Ilmamen |122| . Recently, for the non-time symmetric 
case, Ben-Dov has constructed an example where the apparent horizon does not satisfy 
this inequality |47| . This is not a contradiction with the original Penrose inequality 
which referred to the area of cross-sections of the event horizon, however it does show 
that extending the results beyond time-symmetry would be quite non-trivial. The 'flows' 
which led to the area law in Section O and balance equations in Section [4.2.21 mav be 
potentially useful for this purpose. This approach could lead to an inequality relating 
the area of certain marginally trapped surfaces (the ones connected to future time like 
infinity via a dynamical horizon) to the future limit of the Bondi mass. The framework 
also suggests a program which could shed much light on what John Wheeler called 'the 
issue of the final state': what are the final equilibrium states of a dynamical black hole 
and how, in detail, is this equilibrium reached? From a mathematical perspective an 
important step in addressing this issue is to analyze the non-linear stability of Kerr 
black holes. Consider, then, a neighborhood of the initial data of the Kerr solution in 
an appropriate Sobolev space. One would expect the space-time resulting from evolu- 
tion of this data to admit a dynamical horizon which, in the distant future, tends to an 
isolated horizon with geometry that is isomorphic to that of a Kerr horizon. Can one 
establish that this is what happens? Can one estimate the 'rate' with which the Kerr 
geometry is approached in the asymptotic future? One avenue is to first establish that 
the solution admits a Kerr-Schild type foliation, each leaf of which admits an apparent 
horizon, then show that the world tube of these apparent horizons is a DH, and finally 
study the decay rates of fields along this DH. 

Angular momentum 

As we saw in Section 13 most of the current work on IHs and DHs assumes the presence 
of an axial symmetry 0" on the horizon. A natural question arises: Can one weaken this 
requirement to incorporate situations in which there is only an approximate - rather 
than an exact - symmetry of the horizon geometry? The answer is in the affirmative 
in the following sense. Recall first that the Newman-Penrose component '1/2 is gauge 
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invariant on IHs and DHs. Let F = |^'2p- (While any geometric field could be used here, 
F is the most natural candidate because, for IHs, 4*2 encodes the horizon geometry.) 
If the horizon geometry admits a symmetry, the orbits of the symmetry field (j)"' are 
the level surfaces of F. More generally, let us suppose that the level surfaces of F 
provide a foliation of each good cut S (minus two points) of the horizon. Then, using 
the procedure of section 2.1 of one can introduce on a vector field tangential 
to the foliation, which has the property that it agrees with the symmetry vector field 0" 
whenever the horizon geometry admits a symmetry. The procedure fails if the metric 
on S is at of a round sphere but should work generically. One can then use <j)Q to define 
angular momentum. On IHs this angular momentum is conserved; on DHs it satisfies 
the balance law of Section 14.2.21 and in terms of the initial data, angular momentum 
jiM jias ^jie familiar form (|59|l . IS this proposal viable for numerical simulations? In 
cases ready analyzed, it would be interesting to construct 0q d compare it with the 
symmetry vector field obtained via Killing transport. A better test would be provided 
by non-axisymmetric Brill waves. A second issue associated with angular momentum 
is ether the Kerr inequality J < GIv'P can be violated in the early stages of black hole 
formation or merger, particularly in a non-axisymmetric context. Equations that must 
hold on DHs provide no obvious obstruction 'HTI. Note that such an occurrence is not 
incompatible with the DH finally settling down to a Kerr horizon. For, there is likely 
to be radiation trapped between the DH and the 'peak' of the effective gravitational 
potential that could fall into the DH as time elapses, reducing its angular momentum 
and increasing its mass. The issue of whether a black hole can violate the Kerr inequality 
when it is first formed is of considerable interest to astrophysics j89) . Again, numerical 
simulations involving, say. Brill waves would shed considerable light on this possibility. 

Black hole thermodynamics 

The fact that an integral version of the first law is valid even for non-equilibrium pro- 
cesses, during which the horizon makes a transition from a given state to one which is 
far removed, has interesting thermodynamic ramifications. In non-equilibrium thermo- 
dynamical processes, in general the system does not have time to come to equilibrium, 
whence there is no canonical notion of its temperature. Therefore, while one can still 
interpret the difference E2 — Ei — (work) as the heat AQ absorbed by the system, in 
general there is no longer a clean split AQ — TAS of this term into a temperature 
part and a change in entropy part. If the process is such that the system remains close 
to equilibrium throughout the process, i.e., can be thought of as making continuous 
transitions between a series of equilibrium states, then the difference can be expressed 
as J TdS, where the temperature T varies slowly during the transition. The situation 
on dynamical horizons is analogous. It is only when the horizon geometry is changing 
slowly that the effective surface gravity R of Section [4. 2. 21 would be a good measure of 
temperature, and the horizon area a good measure of entropy (see Section 5.3 of PT]). 
These restricted situations are nonetheless very interesting. Can the black hole entropy 
derivations based on counting of micro-states, such as those of jTJ]], be extended to such 
DHs? In the case of event horizons one would not expect such a procedure to be mean- 
ingful because, as we saw in Eauation l2.2.2l an event horizon can be formed and grow 
in a flat space region in anticipation of a future gravitational collapse. It is difficult to 
imagine how a quasi-local counting of micro-states can account for this phenomenon. 

Perhaps the most surprising aspect of the current status of the theory of black holes is that 
so little is known about their properties in the fully dynamical and non-linear regime of general 
relativity. Indeed, we do not even have a fully satisfactory definition of a dynamical black hole. 
Traditionally, one uses event horizons. But as we discussed in detail, they have several undesirable 
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features. First, they are defined only in space-times where one can unambiguously identify infinity. 
Even in these restricted contexts, event horizons are teleological, can form in a flat region of space- 
time and grow even though there is no flux of energy of any kind across them. When astronomers 
tell us that there is a black hole in the center of Milky Way, they are certainly not referring to 
event horizons. 

Numerical simulations |48l ll()2| suggest that the outermost marginally trapped world-tubes 
become dynamical horizons soon after they are formed. As we saw, dynamical horizons have a 
number of attractive properties that overcome the limitations of event horizons: they are defined 
quasi-locally, can not be formed in flat space-time, and their growth is dictated by balance laws 
with direct physical interpretation. Physically, then, dynamical horizons satisfying the additional 
physical condition £„6(£) < (i.e., SFOTHs) appear to be good candidates to represent the surface 
of a black hole. But so far, our understanding of their uniqueness is rather limited. If a canonical 
dynamical horizon could be singled out by imposing physically reasonable conditions, one could 
use it as the physical representation of an evolving black hole. 

A plausibility argument for the existence of a canonical dynamical horizon was given by Hay- 
ward. Note flrst that on physical grounds it seems natural to associate a black hole with a con- 
nected, trapped region T in space-time (see Section r2.2.1|l . Hay ward |117) sketched a proof that, 
under seemingly natural but technically strong conditions, the dynamical portion of its boundary, 
dT , would be a dynamical horizon H. This H could serve as the canonical representation of the 
surface of an evolving black hole. However, it is not clear whether Hayward's assumptions are 
not too strong. To illustrate the concern, let us consider a single black hole. Then, Hayward's 
argument implies that there are no trapped surfaces outside H. On the other hand there has been 
a general expectation in the community that, given any point in the interior of the event horizon, 
there passes a (marginally) trapped surface through it (see, e.g., This would imply that 

the boundary of the trapped region is the event horizon which, being null, can not qualify as a 
dynamical horizon. However, to our knowledge, this result has not been firmly established. But it 
is clear that this expectation contradicts the conclusion based on Hayward's arguments. Which of 
these two expectations is correct? It is surprising that such a basic issue is still unresolved. The 
primary reason is that very little is known about trapped and marginally trapped surfaces which 
fail to be spherical symmetric. Because of this, we do not know the boundary of the trapped region 
even in the Vaidya solution. 

If it should turn out that the second expectation is correct, one would conclude that Hayward's 
assumptions on the properties of the boundary of the trapped region are not met in physically 
interesting situations. However, this would not rule out the possibility of singling out a canonical 
dynamical horizon through some other conditions (as, e.g., in the Vaidya solution). But since this 
dynamical horizon would not be the boundary of the trapped region, one would be led to conclude 
that in the dynamical and fully non-linear regime, one has to give up the idea that there is a single 
3-manifold that can be interpreted as the black hole surface without further qualifications. For 
certain questions and in certain situations, the dynamical horizon may be the appropriate concept, 
while for other questions and in different situations, the boundary of the trapped region (which 
may be the event horizon) may be more appropriate. 
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